GEOMETRICAL STRUCTURE OF LAPLACIAN EIGENFUNCTIONS 
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Abstract. We summarize the properties of eigenvalues and eigenfunctions of the Laplace opera- 
tor in bounded Euclidean domains with Dirichlet, Neumann or Robin boundary condition. We keep 
the presentation at a level accessible to scientists from various disciplines ranging from mathematics 
to physics and computer sciences. The main focus is put onto multiple intricate relations between 
the shape of a domain and the geometrical structure of eigenfunctions. 
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1. Introduction. This review focuses on the classical eigenvalue problem for the 
Laplace operator A = d'^ /dx\ + ... + ldx\ in an open bounded connected domain 
C M'^ = 2,3,... being the space dimension), 

- Aii^(x) = A^^i^(x) (xG^^), (1.1) 

with Dirichlet, Neumann or Robin boundary condition on a piecewise smooth bound- 
ary dO.'. 



Um{x) = {x e dQ) (Neumann), 2) 



Um{x) = {x e dQ) (Dirichlet) 

d_ 
dn 
d 

^Um(x) + hum(x) = (x G dft) (Robin), 
on 

where d/dn is the normal derivative pointed outwards the domain, and is a posi- 
tive constant. The spectrum of the Laplace operator is known to be discrete, the 
eigenvalues Am are nonnegative and ordered in an ascending order by the index m = 
1,2,3,..., 

(0<)Ai<A2<A3<.../oo (1.3) 

(with possible multiplicities), while the eigenfunctions {um{x)} form a complete basis 
in the functional space L2{^) of measurable and square-integrable functions on Q 
\13S\ I42l| . By definition, the function satisfying Eqs. ( |1.1[ |1.2[ ) is excluded from 
the set of eigenfunctions. Since the eigenfunctions are defined up to a multiplicative 
factor, it is sometimes convenient to normalize them to get the unit L2-norm: 

1/2 

I I dx\Um{x)\A =1 (1.4) 
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(note that there is stih ambiguity up to the multiphcation by e*", with a G M). 

Laplacian eigenfunctions appear as vibration modes in acoustics, as electron wave 
functions in quantum waveguides, as natural basis for constructing heat kernels in 
the theory of diffusion, etc. For instance, vibration modes of a thin membrane (a 
drum) with a fixed boundary are given by Dirichlet Laplacian eigenfunctions 
with the drum frequencies proportional to V%n (418j . A particular eigenmode can 
be excited at the corresponding frequency [44214444] . In the theory of diffusion, an 
interpretation of eigenfunctions is less explicit. The first eigenfunction represents the 
long-time asymptotic spatial distribution of particles diffusing in a bounded domain 
(see below). A conjectural probabilistic representation of higher-order eigenfunctions 
through a Fleming- Viot type model was developed by Burdzy et al. |1QH 1103] . 

The eigenvalue problem |1.2[ ) is archetypical in the theory of elliptic oper- 
ators, while the properties of the underlying eigenfunctions have been thoroughly 
investigated in various mathematical and physical disciplines, including spectral the- 
ory, probability and stochastic processes, dynamical systems and quantum billiards, 
condensed matter physics and quantum mechanics, theory of acoustical, optical and 
quantum waveguides, computer sciences, etc. Many books and reviews were dedicated 
to different aspects of Laplacian eigenvalues, eigenfunctions and their applications 
(see, e.g., [Ill EZl El [Ml [Ml 1123 IWl [1^ The diver- 

sity of notions and methods developed by mathematicians, physicists and computer 
scientists often makes the progress in one discipline almost unknown or hardly acces- 
sible to scientists from the other disciplines. One of the goals of the review is to bring 
together various facts about Laplacian eigenvalues and eigenfunctions and to present 
them at a level accessible to scientists from various disciplines. For this purpose, many 
technical details and generalities are omitted in favor to simple illustrations. While 
the presentation is focused on the Laplace operator in bounded Euclidean domains 
with piecewise smooth boundaries, a number of extensions are relatively straightfor- 
ward. For instance, the Laplace operator can be extended to a second order elliptic 
operator with appropriate coefficients, the piecewise smoothness of a boundary can 
often be relaxed |219[ 1329] , while Euclidean domains can be replaced by Riemannian 
manifolds or weighted graphs [254 . The main emphasis is put onto the geometrical 
structure of Laplacian eigenfunctions and on their relation to the shape of a domain. 
Although the bibliography counts more than five hundred citations, it is far from be- 
ing complete, and readers are invited to refer to other reviews and books for further 
details and references. 

The review is organized as follows. We start by recalling in Sec. [2] general prop- 
erties of the Laplace operator. Explicit representations of eigenvalues and eigenfunc- 
tions in simple domains are summarized in Sec. |3] In Sec. |4]we review the properties 
of eigenvalues and their relation to the shape of a domain: Weyl's asymptotic law, 
isoperimetric inequalities and the related shape optimization problems, and Kac's 
inverse spectral problem. Although eigenfunctions are not involved at this step, valu- 
able information can be learned about the domain from the eigenvalues alone. The 
next step consists in the analysis of nodal lines/surfaces or nodal domains in Sec. [s] 
The nodal lines tell us how the zeros of eigenfunctions are spatially distributed, while 
their amplitudes are still ignored. In Sec. [6| several estimates for the amplitudes of 
eigenfunctions are summarized. Most of these results were obtained during the last 
twenty years. 

Section [7| is devoted to the property of eigenfunctions known as localization. 
We start by recalling the notion of localization in quantum mechanics: the strong 
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localization by a potential (Sec. 7.1), Anderson localization (Sec. 7.2) and trapped 



modes in infinite waveguides (Sec. 7.3). In all three cases, the eigenvalue problem 
is different from Eqs. |1.2| ), due to either the presence of a potential, or the 

unboundness of a domain. Nevertheless, these cases are instructive, as similar effects 
may be observed for the eigenvalue problem |1.2| ). In particular, we discuss in 



Sec. |7.4| an exponentially decaying upper bound for the norm of eigenfunctions in 
domains with branches of variable cross- sectional profiles. Section |7.5| reviews the 
properties of low-frequency eigenfunctions in "dumbbell" domains, in which two (or 
many) subdomains are connected by narrow channels. This situation is convenient for 
a rigorous analysis as the width of channels plays the role of a small parameter [440 . 
A number of asymptotic results for eigenvalues and eigenfunctions were derived, for 
Dirichlet, Neumann and Robin boundary conditions. A harder case of irregular or 
fractal domains is discussed in Sec. |7.6[ Here, it is difficult to identify a relevant 
small parameter to get rigorous estimates. In spite of numerous numerical examples 
of localized eigenfunctions (both for Dirichlet and Neumann boundary conditions), 
a comprehensive theory of localization is still missing. Section |7.7| is devoted to 
high-frequency localization and the related scarring problems in quantum billiards. 
We start by illustrating the classical whispering gallery, bouncing ball and focusing 
modes in circular and elliptical domains. We also provide examples for the case 
without localization. A brief overview of quantum billiards is presented. In the last 
Sec. [8) we mention some issues that could not be included into the review, e.g., 
numerical methods for computation of eigenfunctions or their numerous applications. 

2. Basic properties. We start by recalling basic properties of the Laplacian 
eigenvalues and eigenfunctions (see [7T | I138| 1421] or other standard textbooks). 

(i) The eigenfunctions are infinitely differentiable inside the domain Q. For any 
open subset V C ft, the restriction of Um on V cannot be strictly |3QQJ. 

(ii) Multiplying Eq. ( |1.1| ) by li^, integrating over Q and using the Green's formula 
yield 



J dx \VUm\'^ J dx Um^Q^ iiy ||2 ^ hWn IP 



(2.1) 



where V stands for the gradient operator, and we used Robin boundary condition 



(1.2) in the last equality; for Dirichlet or Neumann boundary conditions, the bound- 
ary integral (second term) vanishes. This formula ensures that all eigenvalues are 
nonnegative. 

(iii) Similar expression appears in the variational formulation of the eigenvalue 
problem, known as the minimax principle \TSS\ 



\Vv\\^r^,^,+_h\\vr 
l|2" 



= minmax ^fl^:^^^^, (2.2) 



where the maximum is over all linear combinations of the form 

V ai0i + ... + am(t>m, 

and the minimum is over all choices of m linearly independent continuous and piecewise- 
differentiable functions (j)m (said to be in the Sobolev space H^{Q)) 1138(1238] . 

Note that the minimum is reached exactly on the eigenfunction Um- For Dirichlet 
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Fig. 2.1. A counter-example for the property of domain monotonicity for Neumann boundary 
condition. Although a smaller rectangle Qi is inscribed into a larger rectangle fl2 (i.e., Qi C ^2), 
the second eigenvalue A2(^^i) = tt^/c^ is smaller than the second eigenvalue X2{^2) = 7r^/a^ (if 
a> b) when c = ^ {a - + (6 - (3^ > a (courtesy by N. Saito; see also j238f . Sec. 1.3.2). 



eigenvalue problem, there is a supplementary condition v = on the boundary dVt 



so that the second term in Eq. (2.2) is canceled. For Neumann eigenvalue problem, 
h = and the second term vanishes again. 

(iv) The minimax principle implies the monotonous increase of the eigenvalues 
with /i, namely \i h < h' ^ then Xm{h) < \m{h'). In particular, any eigenvalue 

Xm{h) of the Robin problem lies between the corresponding Neumann and Dirichlet 
eigenvalues. 

(v) For Dirichlet boundary condition, the minimax principle implies the property 
of domain monotonicity: eigenvalues monotonously decrease when the domain en- 
larges, i.e., Xmi^i) > Xm{^2) if ^1 C 1^2- This property does not hold for Neumann 



or Robin boundary conditions, as illustrated by a simple counter-example on Fig. |2.1 

(vi) The eigenvalues are invariant under translations and rotations of the domain. 
This is a key property for an efficient image recognition and analysis |424[ 1436) 1437] . 
When a domain is expanded by factor a, all the eigenvalues are rescaled by 

(vii) The first eigenfunction ui does not change the sign and can be chosen posi- 
tive. Because of the orthogonality of eigenfunctions, ui is in fact the only eigenfunction 
not changing its sign. 

(viii) The first eigenvalue Ai is simple and strictly positive for Dirichlet and Robin 
boundary conditions; for Neumann boundary condition, Ai = and ui is a constant. 

(ix) The completeness of eigenfunctions in L2{^) can be expressed as 

^Um{x)u^^{y) = 5{x-y) {x.yeVt), (2.3) 

where asterisk denotes the complex conjugate, 5{x) is the Dirac distribution, and the 
eigenfunctions are L2-normalized. Multiplying this relation by a function / G i^2(^) 
and integrating over Vt yields the decomposition of f{x) over Um{x): 



fix) = ^Um{x) / dy f{y) ul^{y). 
h 



(x) The Green function G(x, y) for the Laplace operator which satisfies 

-/^G{x,y) = 8{x-y) (x.yen) (2.4) 

(with an appropriate boundary condition), admits the spectral decomposition over 
the I/2-normalized eigenfunctions 

G{x,y) = Y,X-'u^{x)u*^{y). (2.5) 
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(for Neumann boundary condition, Ai = has to be excluded; in that case, the Green 
function is defined up to an additive constant). 

Similarly, the heat kernel (or diffusion propagator) Gt{x^y) satisfying 

^Gt(x, y) - AGt{x, y) = {x,y e Q), 
Gt=o{x,y) = S{x-y) 
(with an appropriate boundary condition), admits the spectral decomposition 

Gt{x, y) = J2 e-^-'u^{x)ul{y). (2.7) 

m 

The Green function and heat kernel allow one to solve the standard boundary value 
and Cauchy problems for the Laplace and heat equations, respectively |118[ 1139] . 



The decompositions (2.5, 2.7) are the major tool for getting explicit solutions in 
simple domains for which the eigenfunctions are known explicitly (see Sec. |3|. This 
representation is also important for the theory of diffusion due to the probabilistic 
interpretation of Gt(x, y)dx as the conditional probability for Brownian motion started 
at y to arrive in the dx vicinity of x after a time t [SH El lIlZl EOa |249l HQI l42Ql 
1500] . Setting Dirichlet, Neumann or Robin boundary conditions, one can respectively 
describe perfect absorptions, perfect refiections and partial absorption/refiection on 
the boundary [212 . 

For Dirichlet boundary condition, if 1] C Q\ then < Gf^\x^y) < g[^ \x^y) 
|49Qj . In particular, taking Vt' = R^, one gets 



< Gt{x,y) < (4^t)-^/^exp ^-^^J ^ (2-8) 

where the Gaussian heat kernel for free space is written on the right-hand side. The 
above domain monotonicity for heat kernels may not hold for Neumann boundary 
condition [53 . 

(xi) For Dirichlet boundary condition, the eigenvalues vary continuously under 
a "continuous" perturbation of the domain [138J. For Neumann boundary condi- 
tion, the situation is much more delicate. The continuity still holds when a bounded 
domain with a smooth boundary is deformed by a "continuously differentiable trans- 
formation", while in general this statement is false, with an explicit counter-example 
provided in [138 . Note that the continuity of the spectrum is important for numerical 
computations of the eigenvalues by finite element or other methods, in which an irreg- 
ular boundary is replaced by a suitable polygonal or piecewise smooth approximation. 
The underlying assumption that the eigenvalues are very little affected by such do- 
main perturbations, holds in great generality for Dirichlet boundary condition, but is 
much less evident for Neumann boundary condition |lQ5j . The spectral stability of 
elliptic operators under domain perturbations has been thoroughly investigated [105- 
I1Q9[ I226[ 1240] . It is also worth stressing that the spectrum of the Laplace operator in 
a bounded domain with Neumann boundary condition on an irregular boundary may 
not be discrete, with explicit counter-examples provided in [237 . 

3. Eigenbasis for simple domains. We list the examples of "simple" domains, 
in which symmetries allow for variable separations and thus explicit representations 
of eigenfunctions in terms of elementary or special functions. 
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3.1. Intervals, rectangles, parallelepipeds. For rectangle-like domains Q = 
[0,^1] X ... X [O^id] C (with the sizes £i > 0), the natural variable separation yields 

Un,,...,nAxi,...,Xa) = uW(a^l) • ■■ulf}iXd), A„„...,„, = + ... + X^f}, (3.1) 

where the multiple index ni...nd is used instead of m, and uil} (xi) and Xn} {i = 1, (i) 
correspond to the one-dimensional problem on the interval [0,^^]. Depending on 
the boundary condition, Un\x) are sines (Dirichlet), cosines (Neumann) or their 
combinations (Robin): 

uli\x) = sin(7r(n + l)x/£i), X^^^ = 7r^(n + 1)VA^ (Dirichlet), 

u^\x) = cos{7rnx/£i), A^*^ = 7r^n^/-£^, (Neumann), (^3^2) 

u^\x) = sin{anx/£i) + ^ cos(an^/'^i), A^*^ = OL^/£^i, (Robin), 

where n = 0, 1, 2, ... and the coefficients an depend on the parameter h and satisfy the 
equation obtained by imposing the Robin boundary condition in Eq. (1.2) at x = if. 

2ar,. ( ^ 



COS an ^ K^J ^^^^^ ^ ^' (^-^^ 

According to the property (iv) of Sec. [2j this equation has the unique solution an 
on each interval [nyr, (n + l)7r] (n = 0, 1, 2, ...), that makes its numerical computation 
by bisection (or other) method easy and fast. All the eigenvalues Xn^ are simple (not 
degenerate), while 

II (0/ Ml (al^2M,^K'eX^^ , 
\\<\^)\\lm^A)) = i^-^ ^2 -) ' (3-4) 

In turn, the eigenvalues Ani,...,nd can be degenerate if there exists a rational 
ratio (^i/^j)^ (with i 7^ j). For instance, the first Dirichlet eigenvalues of the unit 
square are 27r^, Stt^, Stt^, Stt^, .... , with the twice degenerate second eigenvalue. 
An eigenfunction associated to a degenerate eigenvalue is a linear combination of the 
corresponding functions. For the above example u{xi^X2) = ci sin(7rxi) sin(27rx2) + 
C2 sin(27rxi) sin(7rx2) with any ci and C2 such that + C2 7^ 0. 

3.2. Disk, sector and circular annulus. The rotation symmetry of a circular 
annulus, 1] = {x G : Rq < \x\ < R}^ allows one to write the Laplace operator in 
polar coordinates, 

fxi=rcos(/?, ^^^^1^^^^ 
]x2=rsm(f, dr'^ r dr r'^ d^p'^' 

that leads to variable separation and an explicit representation of eigenfunctions 

Unkl{r,(p) = [Jn{(^nkr/R) + Cn/c^^n (o^n/c^/^)] X \ ^ ^' (3^5) 

\sm{n(f), / = 2 (n 7^ 0), 

where Jn{z) and Yn{z) are the Bessel functions of the first and second kind [U |86l [498] , 
and the coefficients ank and Cnk are set by the boundary conditions ai r = R and 



7 



i?0* 



O^nk 

R 



C^nk 

R 



J'ni^nk) - 



Ro 



h 



R 



) + CnkYni^nk 



Jn{(^nk 

Ro 



R 



Cnk^n{(^nk) 

Ro \ , 



R 



(3.7) 

For each n = 0,1,2,..., the system of these equations has infinitely many solutions 
ank which are enumerated by the index k = 1,2,3,... [498 . The eigenfunctions are 
enumerated by the triple index nkl, with n = 0, 1, 2, ... counting the order of Bessel 
functions, k = 1, 2, 3, ... counting solutions of Eqs. (3.7), and / = 1,2. Since uok2{^^ ^) 
are trivially zero (as sin(n(/p) = for n = 0), they are excluded. The eigenvalues 
\nk = <^n/e/^^' which are independent of the last index are simple for n = and 
twice degenerate for n > 0. In the latter case, an eigenfunction is any nontrivial linear 
combination of Unki and Unk2- The squared L2-norm of the eigenfunction is 

^n,o)R^ \ f 2 



\\Unkl{r,Lp)\\ 



7r(2- 



2at 



2 d2 



h'R 



i^lk 



<k{R) 



(3.8 



where v„fc(r) = Jnia^kr/R) + c„fcy„(a„fer/i?). 

For the special case of a disk (i?o = 0), all the coefRcients Cnk in front of the 
Bessel functions Yn{z) (divergent at 0) are set to 0: 



Unki{r,ip) = Jn{a„kr/R) x 



sin{Tup), 



■ 1, 

2 (n ^ 0), 



(3.9) 



where ank are either the positive roots jnk of the Bessel function Jn{z) (Dirichlet), or 
the positive roots jnk of its derivative Jn{z) (Neumann), or the positive roots of their 
linear combination J!^{z)-\-hJn{z) (Robin). The asymptotic behavior of zeros of Bessel 
functions was thoroughly investigated. For fixed k and large n, the Olver's expansion 
holds jnk - n^6kn^^^^0{n-^^^) (with known coefficients 4) [TMlElIllSTSj, while for 
fixed n and large /c, the McMahon's expansion holds: jnk — 7r(A; + n/2 — 1/4) +0(/c~^) 
[498]. Similar asymptotic relations are applicable for Neumann and Robin boundary 
conditions. 

For a circular sector of radius R and of angle 7r/3, the eigenfunctions are 



Unk{r,(p) = Jn/(3{0Lnkr/R) X 



ir <R, < (/p < 7r/3) 



sin(n(/?//3) (Dirichlet) 
cos{n(p/l3) (Neumann) 

(3.10) 

i.e., they are expressed in terms of Bessel functions of fractional order, and ank are 
the positive roots of Jn/isi^) (Dirichlet) or J'n/jsi^) (Neumann). The Robin boundary 
condition and a sector of a circular annulus can be treated similarly. 

3.3. Sphere and spherical shell. The rotation symmetry of a spherical shell 
in three dimensions, ^7 = {x G : Rq <\x\ < R}^ allows one to write the Laplace 
operator in spherical coordinates. 



xi = r sme/ cos 
X2 = r sin ^ sin 
xs = r cos 6>, 



Qy,2 



2d_ 
r dr 



1 



1 d 



2 V sin (9 9(9 



80 



(3.11) 
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that leads to the variable separation and an explicit representation of eigenfunctions 
Unki{r,0,ip) = [jn{o^nkr/R) ^ Cnkyn{o^nkr/R)]P^{cosO)e'^^, (3.12) 
where jn{z) and yn{z) are the spherical Bessel functions of the first and second kind, 

jn{z) = y^Jn+l/2(^), yn{z) = ]J^Y^^i/2{z) , (3.13) 

P^(z) are associated Legendre polynomials (note that the angular part, P^(cos 6>)e*^'^, 
is also called spherical harmonic and denoted as Yni{0^(p)^ up to a normalization 
factor). The coefficients ank and Cnk are set by the boundary conditions at r = and 



r = Rq similar to Eq. (3.7). The eigenfunctions are enumerated by the triple index 



nkl, with n = 0, 1, 2, ... counting the order of spherical Bessel functions. A: = 1, 2, 3, ... 
counting zeros, and / = — n, —n + 1, ...,n. The eigenvalues Xnk = <^n/c/^^' which are 
independent of the last index /, have the degeneracy 2n + 1. The squared L2-norm of 
the eigenfunction is 



27rR^ 



h^R^ -hR-n{n^l) ]vl^{R) 

(3.14) 



(2n + l)a^, 

al^{Ro/Rf + h^Rl - hRo - n{n + l)Ro/ R]vI^{Ro) 



\\unki{r,0,(p)\\l 



where Vnk{r) = jn{c^nkr/R) + {ankr/R). 

In the special case of a sphere {Rq = 0), one has Cnk = and the equations are 
simplified. For balls and spherical shells in higher dimensions {d > 3), the radial 
dependence of eigenfunctions is expressed through a linear combination of so-called 
ultra-spherical Bessel functions r^~^^'^ Jd_-^_^^{ank'^/ R) and r^~^^'^Yd_-^_^^{ank'^/R)' 

3.4. Ellipse and elliptical annulus. In elliptic coordinates, the Laplace oper- 
ator reads as 

f xi = a cosh r cos l9, ^_ 1 ( ^ ^ ^\ (3 15) 

= asinhrsin6>, (sinh^ r + sin^ l9) V^^^ dO'^ J ' 

where a > is the prescribed distance between the origin and the foci, r > and 
< < 271 are the radial and angular coordinates (Fig. |3.1[ ). An ellipse is a curve 
of constant r = R so that its points (xi,X2) satisfy xf/A^ + ^2/^^ — where R 
is the "radius" of the ellipse and A = a cosh R and B = asinhi? are the major and 
minor semi-axes. Note that the eccentricity e = a/ A = 1/ coshi? is strictly positive. 
A filled ellipse (i.e., the interior of an given ellipse) can be characterized in elliptic 
coordinates as < r < and < < 27r. Similarly, an elliptical annulus (i.e., the 
interior between two ellipses with the same foci) is characterized hy Rq < r < R and 
< 6> < 27r. 

In the elliptic coordinates, the variables can be separated, u{r^9) = f{r)g{0)^ 
from which Eq. ( |1.1[ ) reads as 

jF)d^ + ^cosh(2r)j = - - ^ 

so that both sides are equal to a constant (denoted c). As a consequence, the angular 
and radial parts, g{0) and /(r), are solutions of the Mathieu equation and the modified 
Mathieu equation, respectively [12i 



f{0) + (c - 2(7 cos 20) g{0) = 0, f\r) - {c - 2q cosh 2r) /(r) = 0, 
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B 




A 



Fig. 3.1. Two ellipses of radii R = 0.5 (dashed line) and R = 1 (solid line), with the focal 
distance a = 1. The major and minor semi-axes, A = a cosh R and B = asinhR, are shown by 
black dotted lines. The horizontal thick segment connects the foci. 

where q = and the parameter c is cahed the characteristic value of Mathieu 

functions. Periodic solutions of the Mathieu equation are possible for specific values 
of c. They are denoted as cen(^, q) and sen+i(^, q) (with n = 0, 1, 2, ...) and called the 
angular Mathieu functions of the first and second kind. Each function cq^i (O^q) and 
sen+i(^,^) corresponds to its own characteristic value c (the relation being implicit, 



For the radial part, there are two linearly independent solutions for each char- 
acteristic value c: two modified Mathieu functions Mc^^^ (r, q) and Mc^^^ (r, q) corre- 
spond to the same c as ce^iO^q)^ and two modified Mathieu functions M.s^^}^'^{r^q) 

and y[^\j^i{r.,q) correspond to the same c as sen+i(^,g). As a consequence, there are 
four families of eigenfunctions (distinguished by the index / = 1, 2, 3, 4) in an elliptical 
domain 



where the parameters qnki are determined by the boundary condition. For instance, for 
a filled ellipse of radius R with Dirichlet boundary condition, there are four individual 
equations for the parameter q for each n = 0, 1, 2, ... 



Mc^^\R,qnki) = 0, Mc^^\R,qnk2) = 0, Ms^^l^{R, q^ks) = 0, Ms^^l^{R, q^kA) = 0, 



each of them having infinitely many positive solutions qnki enumerated hy k = 1,2,... 
\T\ 1352] . Finally, the associated eigenvalues are Xnki = Mnki/ct^- 

3.5. Equilateral triangle. Lame discovered the Dirichlet eigenvalues and eigen- 
functions of the equilateral triangle Q = {(xi,X2) G : < xi < 1, < X2 < 
Xi>/3, X2 < a/3(1 — ^i)} by using reflections and the related symmetries |3Qlj : 



see [352] ). 



Unki{r, 0) 

Unk2{r, 0) 

Unks{r, 0) 

Unk4{r, 0) 






(3.16) 
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where 3 divides m -\- m ^ 2n, and n ^ 2m, and the associate eigenfunction is 



Umn{xi,X2) = ^ ± exp 
(m' ,n') 



^(^m'xi + (2n'-m')^)j, (3.17) 



where (m^, n') runs over (— n, m — n), (— n, — m), (n — m, — m), (n — m, n), (m, n) and 
{m^m — n) with the ± sign alternating (see also |339j for basic introduction, as well 
as |315[ I334j ) . Pinsky showed that this set of eigenfunctions is complete in L2{0) 
|395[ [396). Note that the conditions m ^ 2n and n ^ 2m should be satisfied for 
all 6 pairs in the sum that yields one additional condition: m ^ —n. The following 
relations hold: U-m-n = Kn,n^ Un^m = ~Kn,n '^m,o = Um,m' All Symmetric 
eigenfunctions are enumerated by the index (m,0). The eigenvalue corresponds 
to a symmetric eigenfunction if and only if m is a multiple of 3 |395j . 
The eigenfunctions for Neumann boundary condition are 



Umn[xi,X2)= cxp Xi ^ [2n - TU ) J 



{m' ,n') 



(3.18) 



where the only condition is that m-\-n are multiples of 3 (and no sign change). Further 
references and extensions (e.g., to Robin boundary conditions) can be found in a series 
of works by McCartin [342 -345, 347 . McCartin also developed a classification of all 
polygonal domains possessing a complete set of trigonometric eigenfunctions of the 
Laplace operator under either Dirichlet or Neumann boundary conditions [ 346j . 

4. Eigenvalues. 

4.1. Weyl's law. The Weyl's law is one of the first connections between the 
spectral properties of the Laplace operator and the geometrical structure of a bounded 
domain Q. In 1911, Hermann Weyl derived the asymptotic behavior of the Laplacian 
eigenvalues [ 5QH l5Q2j : 

47r^ 

where ^d{^) is the Lebesgue measure of Q. (its area in 2D and volume in 3D), and 

is the volume of the unit ball in d dimensions {T{z) being the Gamma function). As 
a consequence, plotting eigenvalues versus m?^^ allows one to extract the area in 2D 
or the volume in 3D. This result can equivalently be written for the counting function 
A^(A) = #{m : A^ < A} (i.e., the number of eigenvalues smaller than A): 

iV(A) oc ^^1^ X<^'- (A ^ oo). (4.3) 

Weyl also conjectured the second asymptotic term which in two and three dimen- 
sions reads as 

'^AT^^/A {d = 2) 
N(\)o^{ 4;r ^tt (A ^ oo), (4.4) 

M#A3/^T^A (d=3) 
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where /i2(^) and fii{dQ) are the area and perimeter of 1] in 2D, /i3(l^) and ii2{dQ) are 
the volume and surface area of Vt in 3D, and sign "-" (resp. "+") refers to the Dirichlet 
(resp. Neumann) boundary condition. The correction terms which yield information 
about the boundary of the domain, were justified, under certain conditions on Vt (e.g., 
convexity) only in 1980 [2521 [354] (see [11] for a historical review and further details). 
Alternatively, one can study the heat trace (or partition function) 

oo ^ 

Z{t) = I dx Gt{x,x) = e"^"* = / e-^'dN{X) (4.5) 

(here Gt{x,y) is the heat kernel, cf. E g. (|2.6[)), for which the following asymptotic 
expansion holds [SH ESI ETl [103 ESU l^.'^fil 

Z{t) = {AT^t)-"/^ (j^ (^kt'''^ + o(t(^+i)/2)) {t ^ 0), (4.6) 

\fc=0 / 

where the coefRcients Cfe are again related to the geometrical characteristics of the 
domain: 

co = ^id{^), ci = -^fid-i{dn), ... (4.7) 

(see [424 for further discussion). Some estimates for the trace of the Dirichlet Lapla- 
cian were given by Davies |148j (see also |492j for the asymptotic behavior of the heat 
content). 

A number of extensions have been proposed. Berry conjectured that, for irreg- 
ular boundaries, for which the Lebesgue measure in the correction term is infinite, 
the correction term should be A^/^ instead of A^^~^^/^, where H is the Hausdorff 
dimension of the boundary [63l[64j. However, Brossard and Carmona constructed a 
counter-example to this conjecture and suggested a modified version, in which the 
Hausdorff dimension was replaced by Minkowski dimension [88 . The modified Weyl- 
Berry conjecture was discussed at length by Lapidus who proved it for d = 1 [ 303, 304] 
(see these references for further discussion). For dimension d higher than 1, this con- 
jecture was disproved by Lapidus and Pomerance [306 . The correction term to the 
Weyl's formula for domains with rough boundary (in particular, for Lipschitz class) 
was studied by Netrusov and Safarov [368]. Levitin and Vassiliev also considered 
the asymptotic formulas for iterated sets with fractal boundary |320j. Extensions to 
various manifolds and higher order Laplacians were discussed [ 154[ 11551 . 

The high-frequency Weyl's law and the related short-time asymptotics of the heat 
kernel have been thoroughly investigated [11 . The dependence of these asymptotic 
laws on the volume and surface of the domain has found applications in physics. For 
instance, diffusion-weighted nuclear magnetic resonance experiments were proposed 
and conducted to estimate the surface-to- volume ratio of mineral samples and biolog- 
ical tissues [213, 236. 250i MB M [Ml [351 [454, . 

The multiplicity of eigenvalues is yet a more difficult problem [363 . From basic 
properties (see Sec. [2|, the first eigenvalue Ai is simple. Cheng proved that the 
multiplicity m{\2) of the second Dirichlet eigenvalue A2 is not greater than 3 [128^ . 
This inequality is sharp since an example of domain with m{\2) = 3 was constructed. 
For /c > 3, Hoffmann-Ostenhof et al. proved the inequality m{Xk) <2k — 3 [ 24611247] . 
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4.2. Isoperimetric inequalities for eigenvalues. In the low- frequency limit, 
the relation between the shape of a domain and the associated eigenvalues manifests in 
the form of isoperimetric inequalities. Since there are many excellent reviews on this 
topic, we only provide a list of the best-known inequalities, while further discussion 
and references can be found in [271 El ESI El ^221 EMI MD\ SOll \MM • 

(i) The Rayleigh-Faber-Krahn inequality states that the disk minimizes the first 
Dirichlet eigenvalue Ai among all planar domains of the same area /i2(^), i-e. 

where jj^^i is the first positive zero of Ju{z) (e.g., jo,i ~ 2.4048...). This inequality 
was conjectured by Lord Rayleigh and proven independently by Faber and Krahn 
|174[ I289j . The corresponding isoperimetric inequality in d dimensions. 



was proven by Krahn |29Qj . 

Another lower bound for the first Dirichlet eigenvalue for a simply connected 
planar domain was obtained by Makai [333^ and later rediscovered (in a weaker form) 
by Hayman [231] 

Af > 4 (4.10) 



where a is a constant, and 



p = maxmin{|x-^|} (4.11) 



is the inradius of Q (i.e., the radius of the largest ball inscribed in Q). The above 
inequality means that the lowest frequency (bass note) can be made arbitrarily small 
only if the domain includes an arbitrarily large circular drum (i.e., p goes to infinity). 
The constant a, which was equal to 1/4 in the original Makai's proof (see also [380') 
and to 1/900 in the Hayman's proof, was gradually increased, to the best value (up 
to date ) a = 0.6197... by Banuelos and Carroll [39]. For convex domains, the lower 



bound ( |4.10[ ) with a = 7r^/4 « 2.4674 was derived much earlier by Hersch [241j, with 
the equality if and only if Q is an infinite strip (see also a historical overview in [29]). 

An obvious upper bound for the first Dirichlet eigenvalue can be obtained from 
the domain monotonicity (property (v) in Sec. [2|: 

Af < Af (Bp) = (4.12) 

with the first Dirichlet eigenvalue Xf{Bp) for the largest ball Bp inscribed in Q {p is 
the inradius). However, this upper bound is not accurate in general. Polya and Szego 
gave another upper bound for planar star-shaped domains [401 . Freitas and Krejcifik 
extended their result to higher dimensions [192 : for a bounded strictly star-shaped 
domain 1] C with locally Lipschitz boundary, they proved 



(4.13) 
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where the function F{rt) is defined in jl92j . From this inequality, they also deduced 
a weaker but more explicit upper bound which is applicable to any bounded convex 
domain in W^: 

Af < Af (50 f^^. (4.14) 

The second Dirichlet eigenvalue X2 is minimized by the union of two identical 
balls 

This inequality, which can be deduced by looking at nodal domains for U2 and using 






Rayleigh-Faber-Krahn inequality (4.9) on each nodal domain, was first established by 
Krahn |29Qj . It is also sometimes attributed to Peter Szego (see [404 ). Note that 
finding the minimizer of X2 among convex planar sets is still an open problem |239J. 
Bucur and Henrot proved the existence of a minimizer for the third eigenvalue in the 
family of domains in of given volume, although its shape remains unknown [93]. 
The range of the first two eigenvalues was also investigated [92] 1505] . 

The first nontrivial Neumann eigenvalue (as = 0) also satisfies the isoperi- 
metric inequality 

A^<(7^) ih.f, (4.16) 



which states that is maximized by a d-dimensional ball (here jj^^i is the first positive 
zero of the function ■j^[z'^~^^'^ Jd_-^_^^{z)] which reduces to Jl{z) and Y^2/7r jl{z) for 
d = 2 and d = 3, respectively). This inequality was proven for simply-connected 
planar domains by Szego [479 and in higher dimensions by Weinberger [499]. Polya 
conjectured the following upper bound for all Neumann eigenvalues |4Q3j in planar 
bounded regular domains (see also [450]) 

A:<fc^ (n = 2,3,4,...) (4.17) 

(the domain is called regular if its Neumann eigenspectrum is discrete, see [204 for 
details). This inequality is true for all domains that tile the plane, e.g., for any triangle 



and any quadrilateral [405 . For n = 2, the inequality (4.17) follows from (4.16). For 
n > 3, Polya's conjecture is still open, although Kroger proved a weaker estimate 
A^ < 87r(n — 1) [292]. Recently, Girouard et al. obtained a sharp upper bound for 
the second nontrivial Neumann eigenvalue A3^ for a regular simply-connected planar 
domain [204j : 

A? , ,4.1s, 

with the equality attained in the limit by a family of domains degenerating to a 
disjoint union of two identical disks. 

Payne and Weinberger obtained the lower bound for the second Neumann eigen- 
value in d dimensions [389^ 



TT 



2 



A^^ > ^, (4.19) 
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where S is the diameter of 

S = max {\x - y\}. (4.20) 

This is the best bound that can be given in terms of the diameter alone in the sense 
that (5^ tends to tt^ for a parahelepiped ah but one of whose dimensions shrink to 



Szego and Weinberger noticed that Szego's proof of the inequahty (4.16) for planar 
simply connected domains extends to prove the bound 

J_ ^ (421) 

with equality if and only if ^7 is a disk |479[ I499j . Ashbaugh and Benguria derived 
another bound for arbitrary bounded domain in [25] 

In particular, one gets I/X2 + 1/^3^ ^ for d = 2 (see also extensions in j244[ 

[506]). 

(ii) The Payne- Poly a- Weinberger inequality^ which can also be called Ashbaugh- 
Benguria inequality^ concerns the ratio between first two Dirichlet eigenvalues and 
states that 

with equality if and only if is the d-dimensional ball. This inequality (in 2D form) 
was conjectured by Payne, Polya and Weinberger [388 and proved by Ashbaugh and 
Benguria in 1990 [23ti26]. A weaker estimate /Af < 1 + 4/^/ was proved for d = 2 



in the original paper by Payne, Polya and Weinberger j388j . 

(iii) Singer et al. derived the upper and lower estimates for the spectral (or 
fundamental) gap between the first two Dirichlet eigenvalues for a smooth convex 
bounded domain Q in (in [465], a more general problem in the presence of a 
potential was considered): 

^>A--Af>^, (4.24) 



where S is the diameter of Q and p is the inradius |465j . For a convex planar domain, 
Donnelly proposed a sharper lower estimate [1? 



Af-Af>^. (4.25) 

However, Ashbaugh et al. pointed out on a flaw in the proof [30^. The estimate was 
later rigorously proved by Andrews and Clutterbuck for any bounded convex domain 
in R'^, even in the presence of a semi-convex potential [10 (for more background 
on the spectral gap, see notes by Ashbaugh |28)). 
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(iv) The isoperimetric inequalities for Robin eigenvalues are less known. Daners 
proved that among all bounded domains Q C of the same volume, the ball B 
minimizes the first Robin eigenvalue [95l I144j 



X^{n) > Af (5). (4.26) 

Kennedy showed that among all bounded domains in R^, a domain B2 composed of 
two disjoint balls minimizes the second Robin eigenvalue [278j 



X^{n) > (^2). (4.27) 

(v) The minimax principle ensures that the Neumann eigenvalues are always 
smaller than the corresponding Dirichlet eigenvalues: A^ < A^. Polya proved A2^ < 
Af |4Q2j while Szego got a sharper inequality A2^ < cXf for a planar domain bounded 
by an analytic curve, where c = (ji,i/jo,i)^ ^ 0.5862... [479 (note that this result also 



follows from inequalities (4.8, 4.16)). Payne derived a stronger inequality for a planar 
domain with a boundary- A^^2 < ^ [387 - Levine and Weinberger 

generalized this result for higher dimensions d and proved that A^^^ < A^ for all 
n when ft is smooth and convex, and that A^^^ < if Q is merely convex [319 . 
Friedlander proved the inequality A^^^ < A^ for a general bounded domain with a 
boundary |194j . Filonov found a simpler proof of this inequality in a more general 
situation (see |187j for details). 

Many other inequalities can be found in several reviews [271 EHl EH] . It is worth 
noting that isoperimetric inequalities are related to shape optimization problems ^ 
[851 IMl [TT2l[TTm [3971 [468] , 

4.3. Kac's inverse spectral problem. The problem of finding relations be- 
tween the Laplacian eigenspectrum and the shape of a domain was formulated in the 
famous Kac's question "Can one hear the shape of a drum?" [266 . In fact, the drum's 
frequencies are uniquely determined by the eigenvalues of the Laplace operator in the 
domain of drum's shape. By definition, the shape of the domain fully determines the 
Laplacian eigenspectrum. Is the opposite true, i.e., does the set of eigenvalues which 
appear as "fingerprints" of the shape, uniquely identify the domain? The negative an- 
swer to this question for general planar domains was given by Gordon and co-authors 
[20T who constructed two different (nonisometric) planar polygons (Fig. |4.1^ ,b) with 
the identical Laplacian eigenspectra, both for Dirichlet and Neumann boundary con- 
ditions (see also [59]). Their construction was based on Sunada's paper on isospectral 
manifolds [478]. An elementary proof, as well as many other examples of isospectral 



domains, were provided by Buser and co-workers |lllj and by Chapman [124 (see 



Fig. 4.1 3, d). An experimental evidence for this not "hearing the shape" of drums 
was brought by Sridhar and Kudrolli [471^ (see also [135 ). In all these examples, 
isospectral domains are either non-convex, or disjoint. Gordon and Webb addressed 
the question of existence of isospectral convex connected domains and answered this 
question positively (i.e., negatively to the original Kac's question) for domains in Eu- 
clidean spaces of dimension d > 4 |2Q8j . To our knowledge, this question remains 
open for convex domains in two and three dimensions, as well as for domains with 
smooth boundaries. It is worth noting that the positive answer to Kac's question 
can be given for some classes of domains. For instance, Zelditch proved that for do- 
mains that possess the symmetry of an ellipse and satisfy some generic conditions on 
the boundary, the spectrum of the Dirichlet Laplacian uniquely determines the shape 
|51Qj . Later, he extended this result to real analytic planar domains with only one 
symmetry [5 111 1512j . 
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Fig. 4.1. Two examples of nonisometric domains with the identical Laplace operator eigen- 
spectra (with Dirichlet or Neumann boundary conditions): the original example (shapes ^a^ and 'b') 
constructed by Gordon et al. \207i , and a simpler example with disconnected domains (shapes ^c^ 
and 'd') by Chapman \124i - ^'^ latter case, the eigenspectrum is simply obtained as the union 
of the eigenspectra of two subdomains known explicitly. For instance, the Dirichlet eigenspectrum 
is {7r2(m2 : m, n G N} U {7r2((z/2)2 + (j72)2) : i,j G N, i > j} . 

A somewhat similar problem was recently formulated for domains in which one 
part of the boundary admits Dirichlet boundary condition and the other Neumann 
boundary condition. Does the spectrum of the Laplace operator determine uniquely 
which condition is imposed on which part? Jakobsen and co-workers gave the negative 
answer to this question by assigning Dirichlet and Neumann conditions onto different 
parts of the boundary of the half-disk (and some other domains), in a way to produce 
the same eigenspectra [255 . 

The Kac's inverse spectral problem can also be seen from a different point of 
view. For a given sequence < Ai < A2 < A3 < whether does exist a domain Vt 
in for which the Laplace operator with Dirichlet or Neumann boundary condition 
has the spectrum given by this sequence. A similar problem can be formulated for a 
compact Riemannian manifold with arbitrary Riemannian metrics. Colin de Verdiere 
studied these problems for finite sequences {An}^=i and proved the existence of such 
domains or manifolds under certain restrictions |137j . We also mention the work of 
Sleeman who discusses the relationship between the inverse scattering theory (i.e., the 
Helmholtz equation for an exterior domain) and the Kac's inverse spectral problem 
(i.e., for an interior domain) |466j (see |133j for further discussion on inverse eigenvalue 
problems). 

5. Nodal lines. The first insight onto the geometrical structure of eigenfunc- 
tions can be gained from their nodal lines. Kuttler and Sigillito gave a brief overview 
of the basic properties of nodal lines for Dirichlet eigenfunctions in two dimensions 
|3QQj that we partly reproduce here: 

"The set of points in Vt where Ujn = is the nodal set of Um- By the unique 
continuation property, it consists of curves that are in the interior of Vt. Where 
nodal lines cross, they form equal angles |138j . Also, when nodal lines intersect 
a portion of the boundary, they form equal angles. Thus, a single nodal line 
intersects the boundary at right angles, two intersect it at 60° angles, and so forth. 
Courant's nodal line theorem [138 states that the nodal lines of the m-th eigenfunction 
divide Vt into no more than m subregions (called nodal domains): Vrn ^ being 
the number of nodal domains. In particular, ui has no interior nodes and so Ai is a 
simple eigenvalue (has multiplicity one)." 

It is worth noting that any eigenvalue A^ of the Dirichlet-Laplace operator in Vt is 
the first eigenvalue for each of its nodal domains. This simple observation allows one 
to construct specific domains with a prescribed eigenvalue (see |3QQj for examples). 
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Fig. 5.1. The nodal lines of a Dirichlet eigenfunction u(xi,X2) on the unit square, with the 
associated eigenvalue A = 55257r^ of multiplicity 12. The eigenfunction was obtained as a linear 
combination of terms sin(7mia::i) sin(7rn2a:2); with n\-\-n^ = 5525 and randomly chosen coefficients. 
For comparison, another eigenfunction with the same eigenvalue, sin(507ra;i) sin(557ra;2); is shown. 



Eigenfunctions with few nodal domains were constructed in |138[ 1321] . 

Even for such a simple domain as a square, the nodal lines and domains may 
have complicated structure, especially for high-frequency eigenfunctions (Fig. 5.1). 
This is particularly true for degenerate eigenfunctions for which one can "tune" the 
coefficients of the corresponding linear combination to modify continuously the nodal 
lines. 

Pleijel sharpened the Courant's theorem by showing that the upper bound m for 
the number of nodal domains is attained only for a finite number of eigenfunctions 
[399 . Moreover, he obtained the upper limit: lim Vrn/'^ = 4/jo i ^ 0.691.... Note 
that Lewy constructed spherical harmonics of any degree n whose nodal sets have one 
component for odd n and two components for even n implying that no non-trivial 
lower bound for Vrn is possible |321j . We also mention that the counting of nodal 
domains can be viewed as partitioning of the domain into a fixed number of subdo- 
mains and minimizing an appropriate "energy" of the partition (e.g., the maximum 
of the ground state energies of the subdomains) . When a partition corresponds to an 
eigenfunction, the ground state energies of all the nodal domains are the same, i.e., it 
is an equipartition pO^. 

Blum et al. considered the distribution of the (properly normalized) number 
of nodal domains of the Dirichlet-Laplacian eigenfunctions in 2D quantum billiards 
and showed the existence of the limiting distribution in the high-frequency limit (i.e., 
when oo) [76 . These distributions were argued to be universal for systems with 

integrable or chaotic classical dynamics that allows one to distinguish them and thus 
provides a new criterion for quantum chaos (see Sec. 7.7.4). It was also conjectured 
that the distribution of nodal domains for chaotic systems coincides with that for 
Gaussian random functions. 

Bogomolny and Schmit proposed a percolation-like model to describe the nodal 
domains which permitted to perform analytical calculations and agreed well with nu- 
merical simulations [77] . This model allows one to apply ideas and methods developed 
within the percolation theory |472j to the field of quantum chaos. Using the analogy 
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with Gaussian random functions, Bogomolny and Schmit obtained that the mean and 
variance of the number of nodal domains grow as m, with exphcit formulas for the 
prefactors. From the percolation theory, the distribution of the area s of the connected 
nodal domains was conjectured to follow a power law, n{s) ex ^-i^t/qi^ confirmed 
by simulations [77]. In the particular case of random Gaussian spherical harmonics, 
Nazarov and Sodin rigorously derived the asymptotic behavior for the number of 
nodal domains of the harmonic of degree n \367\. They proved that as n grows to 
infinity, the mean of Vnl^ tends to a positive constant, and that Vnl^ exponentially 
concentrates around this constant (we recall that the associate eigenvalue is n(n + l)). 

The geometrical structure of nodal lines and domains has been intensively stud- 
ied (see |366[ 1400] for further discussion of the asymptotic nodal geometry). For 
instance, the length of the nodal line of an eigenfunction of the Laplace operator in 
two-dimensional Riemannian manifolds was separately investigated by Briining, Yao 
and Nadirashvili who obtained its lower and upper bounds [90[ I362[ l5Q7j . In addi- 
tion, a number of conjectures about the properties of particular eigenfunctions were 
discussed in the literature. We mention three of them: 

(i) In 1967, Payne conjectured that the second Dirichlet eigenfunction U2 cannot 
have a closed nodal line in a bounded planar domain |39Q[ 1394] . This conjecture was 
proved for convex domains O [353 and disproved by non-convex domains [245 , see 
also [2T61I256] . 

(ii) The hot spots conjecture formulated by J. Ranch in 1974 says that the 
maximum of the second Neumann eigenfunction is attained at a boundary point. 
This conjecture was proved by Banuelos and Burdzy for a class of planar domains 
[4T] but in general the statement is wrong, as shown by several counter-examples 



(iii) Liboff formulated several conjectures; one of them states that the nodal 
surface of the first-excited state of a 3D convex domain intersects its boundary in a 
single simple closed curve [323]. 

The analysis of nodal lines that describe zeros of eigenfunctions, can be extended 
to other level sets. For instance, a level set of the first Dirichlet eigenfunction u\ 
on a bounded convex domain 1] G is itself convex |274j . Grieser and Jerison 
estimated the size of the first eigenfunction uniformly for all convex domains |217j . 
In particular, they located the place where u\ achieves its maximum to within a 
distance comparable to the inradius, uniformly for arbitrarily large diameter. Other 
geometrical characteristics (e.g., the volume of a set on which an eigenfunction is 
positive) can also be analyzed [364 . 

6. Estimates for Laplacian eigenfunctions. The "amplitudes" of eigenfunc- 
tions can be characterized either globally by their norms 



or locally by pointwise estimates. Since eigenfunctions are defined up to a multiplica- 
tive constant, one often uses I/2(^) normalization: ||ix||2 = 1. Note also the limiting 
case of Loo-norm 



(the first equality is the general definition, while the second equality is applicable 
for eigenfunctions). It is worth recalling Holder's inequality for any two measurable 



isi una una ESI]. 




(6.1) 




OO — 



esssup^^j^ |'^(^)| = niax \u{x)\ 



(6.2) 
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functions u and v and for any positive p, q such that l/p-\-l/q = 1: 

\\uv\U < \\u\\4v\U. (6.3) 

For a bounded domain ^7 C M'^ (with a finite Lebesgue measure /j.d{^))^ Holder's 
inequality implies 

\\u\\p<[iid{af^-V\\u\\p, il<p<p'). (6.4) 
We also mention Minkowski's inequality for two measurable functions and any p > 1: 

\\u + v\\p<\\u\\p+\\v\\p. (6.5) 

6.1. First (ground) Dirichlet eigenfunction. The Dirichlet eigenfunction ui 
associated with the first eigenvalue Ai > does not change the sign in Q and may be 
taken to be positive. It satisfies the following inequalities. 

(i) Payne and Rayner showed in two dimensions that 

|K||2<^ ll^xilli, (6.6) 
V47r 



with equality if and only if is a disk |39H I392j . Kohler- Jobin gave an extension of 
this inequality to higher dimensions |283j (see |132( I284( 1392] for other extensions): 

hih< , L^ =K||i. (6.7) 

(ii) Payne and Stakgold derived two inequalities for a convex domain in 2D 

^IKI|.<lklU (6.8) 

and 

u^{x) < \x-d^\^^\\m\\r (xGl^), (6.9) 

where \x — d^\ is the distance from a point x in ^7 to the boundary dVt [393". 

(iii) van den Berg proved the following inequality for L2-normalized eigenfunction 
ui when is an open, bounded and connected set in (d = 2, 3, . . . ): 

d-2 

with equality if and only if (7 is a ball, where p is the inradius (Eq. ( 4.11| )) [493^ . 
van den Berg also conjectured the stronger inequality for an open bounded convex 
domain C M^: 

\\ui\\^<Cdp-'"\pl5f'\ (6.11) 
where 5 is the diameter of Vt^ and Cd is a universal constant independent of Vt. 
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(iv) Kroger obtained another upper bound for \\ui \ \qq for a convex domain 1^ C M. . 
Suppose that Xi{D) > A{S) for every convex subdomain D C with fid{D) < Sjj.d{^) 
and positive numbers S and A{5). The first eigenfunction ui which is normalized such 
that WuiWl = satisfies 



1^1 Hoc < CdS-'^^ [1 + In \\u, lU - ln(l - Xi/A{S))] 



d/2 



(6.12) 



with a universal positive constant Cd which depends only on the dimension d |293j . 

(v) Pang investigated how the first Dirichlet eigenvalue and eigenfunction would 
change when the domain slightly shrinks |383[ I384j . For a bounded simply connected 
open set 1] C , let 

n,^{x en:\x-dn\> e} 

be its interior, i.e., Q without an e boundary layer. Let and be the Dirichlet 
eigenvalues and L2-normalized eigenfunctions in ft^ (with A^ = Xm and = Um 
referring to the original domain Q). Then, for all e G (0,p/2), 

|A^i - All <Ci 6^/2, 
- T,u\ lU^(o) < \C2 + C3(A2 - Ai)-i/2 + C4(A2 - Ai)-i' 



.1/2 



(6.13) 



where p is the inradius of (Eq. ( |4.11[ )), is the extension operator from Q,^ to f2, 
and 



C4 



-3/2^9/4 



37r9/4' 
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-5/2^4 
-7/2^7 



C2 
1 - 



3/2013/4 



9j 



2^ 

^15/4 ' 



8IV2 



a n 



15/2 



7r3/4' 

1 + 187i;93/4 , 



8I7? 
7.3/2 



^3/2 



where /3 = /i2(^)/p^, <^ is the constant from Eq. (4.10) (for which one can use the 
best known estimate a = 0.6197... from [39]), and 71 and 72 are the first and second 
Dirichlet eigenvalues for the unit disk: 71 = j^ ^ ^ 5.7832 and 72 = Ji 1 ^ 14.6820. 
Moreover, when Q is the cardioid in R^, the term e^/^ cannot be improved}^ 

In addition, Davies proved that for a bounded simply connected open set 1] G 
and for any /3 G (0, 1/2), there exists > 1 such that jl49j 



lAl-Ail <c^e'' (6.14) 

for all sufficiently small e > 0. The estimate also holds for higher Dirichlet eigenvalues. 
6.2. Estimates applicable for all eigenfunctions. 

6.2.1. Estimates through the Green function. Using the spectral decom- 
position (2.5) of the Green function G{x^y), one can rewrite Eq. ( |1.1| ) as 

{x) = Xm / G{x,y)um{y)dy, 



^ In the original paper |384| . the coefficient C4 in Eq. (1.5) should be multiplied by the omitted 
prefactor v^|^| that follows from the derivation. 
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from which Holder inequality (6.3) yields a family of simple pointwise estimates 

\\^ (j \G{x,yWdy\ , (6.15) 

with any p > 1. Here, a single function of x in the right-hand side bounds all the 
eigenfunctions. In particular, for p = 1, one gets 

(6.16) 



{x)\ < A 

m r^m Woo 

For Dirichlet boundary condition, G(x, y) is positive everywhere in Q so that 

m r^m CO 

U{x), U{x)= [ G{x,y)dy, 



(6.17) 



where U{x) solves the boundary value problem 

- AU{x) = 1 (a; e il), U{x) 



(x e dQ). 



(6.18) 

The solution of this equation is known to be the mean first passage time to the 
boundary dfl from an interior point x [420]. The inequalities (6.16, 6.17[ ) (or their 
extensions) were reported by Moler and Payne [360 (Sect. |6.2.2 ) and were used by 
Filoche and Mayboroda for determining the geometrical structure of eigenfunctions 

[ng (Sect. 



6.2.6). Note that the function U{x) was also considered by Gorelick et 



al. for a reliable extraction of various shape properties of a silhouette, including part 
structure and rough skeleton, local orientation and aspect ratio of different parts, and 
convex and concave sections of the boundaries [209 . 

6.2.2. Bounds for eigenvalues and eigenfunctions of symmetric oper- 
ators. Moler and Payne derived simple bounds for eigenvalues and eigenfunctions 
of symmetric operators by considering their extensions |360j . As a typical example, 
one can think of the Dirichlet-Laplace operator in a bounded domain (symmet- 
ric operator A) and of the Laplace operator without boundary conditions (extension 
A^). An approximation to an eigenvalue and eigenfunction of A can be obtained by 
solving a simpler eigenvalue problem A^u^ = X^u^ without boundary condition. If 
there exists a function w such that A^w = {) and = ix* at the boundary of 17 and if 



ll^*IU2(^7) 



< 1, then there exists an eigenvalue \k of A satisfying 



l+e 



< |Afc| < 



|A*| 



(6.19) 



Moreover, if ||l2(^7) = 1 '^fc is the L2-normalized projection of onto the 
eigenspace of A^, then 

1/2 



^* - Uk\\L2in) 



< 



1 



with 



mm 



I A, 



(6.20) 



If is a good approximation to an eigenfunction of the Dirichlet-Laplace oper- 
ator, then it must be close to zero on the boundary of 17, yielding small e and thus 



accurate lower and upper bounds in (6.19). The accuracy of the bound (6.20) also 



depends on the separation a between eigenvalues. 

In the same work, Moler and Payne also provided pointwise bounds for eigenfunc- 



tions that rely on Green's functions (an extension of Sec. 6.2.1). 
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6.2.3. Estimates for Lp-norms. Chiti extended Payne- Ray ner's inequality 
(6.6) to the eigenfunctions of linear elliptic second order operators in divergent form, 



with Dirichlet boundary condition |132j . For the Laplace operator in a bounded 
domain 17 C M^, Chiti's inequality for any real numbers q>p>0 states: 



ll^ll, < 11^.11, (du;,)i-^AS^^J-^- ^ (6.21) 

^ / 



where ujd is given by Eq. (4.2). 

6.2.4. Pointwise bounds for Dirichlet eigenfunctions. Banuelos derived a 
pointwise upper bound for L2-normalized Dirichlet eigenfunctions [40] 

Mx)\<\T (xen). (6.22) 

van den Berg and Bolthausen proved several estimates for L2 -normalized Dirichlet 
eigenfunctions [491 . Let 1] C = 2, 3, . . . ) be an open bounded domain with 

boundary which satisfies an a-uniform capacitary density condition with some 
a e (0, 1], i.e. 

Csip{dn n B{x; r)} > aCap{5(x; r)}, x e dQ, < r < (6.23) 

where B{x^ r) is the ball of radius r centered at x, S is the diameter of Q (Eq. ( |4.2Q )), 
and Cap is the logarithmic capacity for d = 2 and the Newtonian (or harmonic) 
capacity for d > 2 (the harmonic capacity of an Euclidean domain presents a measure 
of its "size" through the total charge the domain can hold at a given potential energy 
|161j ). The condition ( 6.23| ) guarantees that all points of are regular. The following 
estimates hold 

(i) in two dimensions {d = 2), for all m = 1,2,... and all x G Q such that 
\x — ^Q\^/Xm < 1, one has 



^ln{\x-dQ\VX _^ ^ 
(ii) in higher dimensions {d > 2), for all m = 1,2,... and all x € such that 

/ 6 \ l+7(<i-l)/(d-2) 

\x-dn\y^<(^] , (6.25) 



"^^^^= ln(2(2/a)V(^-^)) '°"^ 



/ X i((l/7) + (d-l)/(d-2))-^ 

\uM\<^>^T(\x-^n\^/X^] . (6.26) 
(iii) for a planar simply connected domain and all m = 1, 2, . . . , 

\um{x)\ < m 29/2^1/4 (M^)^ 1^ _ dn\i/2 ^ (6.27) 
where p is the inradius of Q (see Eq. ( 4.11[ )), and the inequality is sharp. 
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6.2.5. Upper and lower bounds for normal derivatives of Dirichlet eigen- 
functions. Suppose that M is a compact Riemannian manifold with boundary and 
u is an L2-normahzed Dirichlet eigenfunction with eigenvalue A. Let be its normal 
derivative at the boundary. A scaling argument suggests that the I/2-norm of ip will 
grow as a/a as A ^ oo. Hassell and Tao proved that 

Cm ^ IIV^IIlsCOM) ^ C'mV^, (6.28) 

where the upper bound holds for any Riemannian manifold, while the lower bound 
is valid provided that M has no trapped geodesies |228j . The positive constants cm 
and Cm depend on M, but not on A. 

6.2.6. Estimates for restriction onto a subdomain. For a bounded domain 
1] C M^, Filoche and Mayboroda obtained the upper bound for the L2-norm of a 
Dirichlet-Laplacian eigenfunction u associated to A, in any open subset D C Q [186 : 

< (^1 + lbllL,(D), (6.29) 

where the function v solves the boundary value problem in D: 

Av = (xeD), v = u (xedD). 

and is the distance from A to the spectrum of the Dirichlet-Laplace operator 

in D. Note also that the above bound was proved for general self- adjoint elliptic 



operators [186 . When combined with Eq. (6.17), this inequality helps one investigate 



the spatial distribution of eigenfunctions because harmonic functions are in general 
much easier to analyze or estimate than eigenfunctions. 

We complete the above estimate by a lower bound |371j 



\u\\L,iD)>J^^^^\\v\\L,iD), (6.30) 



where Xi{D) is the first Dirichlet eigenvalue of the subdomain D. 

7. Localization of eigenfunctions. "Localization" is defined in the Webster's 
dictionary as "act of localizing, or state of being localized" . The notion of localization 
appears in various fields of science and often has different meanings. Throughout this 
review, a function u defined on a domain 1] C M^, is called Lp-localized (for p > 1) if 
there exists a bounded subset l^o C 1^ which supports almost all L^-norm of i.e. 

"""^-^^^^-^ « 1 and ^«L (7.1) 

Qualitatively, a localized function essentially "lives" on a small subset of the domain 
and takes small values on the remaining part. For instance, a Gaussian function 
exp(— x^) on = R is L^-localized for any p > 1 since one can choose Qq = [—a, a] 
with large enough a so that the ratio of L^-norms can be made arbitrarily small, while 
the ratio of lengths /ii(l]o)/Mi(^) is strictly 0. In turn, when Q = [—A, A], the local- 
ization character of exp(— x^) on Q becomes dependent on A and thus conventional. 



A simple calculation shows that both ratios in (7.1) cannot be simultaneously made 
smaller than 1/{A + 1) for any p > 1. For instance, if A = 3 and the "threshold" 1/4 
is viewed small enough, then we are justified to call exp(— x^) localized on [—A, A]. 
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This example illustrates that the above inequalities do not provide a universal quanti- 
tative criterion to distinguish localized from non-localized (or extended) functions. In 
this section, we will describe various kinds of localization for which some quantitative 
criteria can be formulated. We will also show that the choice of the norm (i.e., p) 
may be important. 

Another "definition" of localization was given by Felix et at. who combined L2 
and I/4 norms to define the "existence area" as jl75j 

A function u was called localized when its existence area S{u) was much smaller than 
the area /i2(^) |175j (this definition trivially extends to other dimensions). In fact, if 
a function is small in a subdomain, the fourth power diminishes it stronger than the 
second power. For instance, if 1] = (0, 1) and is 1 on the subinterval VLq = (1/4, 1/2) 
and otherwise, one has ||i^||l2(^7) = 1/2 and ||ii||L4(Q) = 1/a/2 so that S{u) = 1/4, 
i.e., the length of the subinterval Qq. This definition is still qualitative: e.g., in the 
above example, is the ratio S{u)/ = 1/4 small enough to call u localized? Note 
that a whole family of "existence areas" can be constructed by comparing Lp and Lq 
norms (with p < q)^ 

S„W=({|^)-- (7.3) 

7.1. Bound quantum states in a potential. The notion of bound, trapped or 
localized quantum states is known for a long time [TH 1421] . The simplest "canonical" 
example is the quantum harmonic oscillator, i.e., a particle of mass m in a harmonic 
potential of frequency u which is described by the Hamiltonian 

where p = —ihdx is the momentum operator, and f = x is the position operator {h 
being the Planck's constant). The eigenfunctions of this operator are well known: 

= V ^ l^i \^ ) H-Wr^M^ (7.5) 

where Hn{x) are the Hermite polynomials. All these functions are localized in a 
region around the minimum of the harmonic potential (here, x = 0), and rapidly 



decay outside this region. For this example, the definition (7.1) of localization is 
rigorous. In physical terms, the presence of a strong potential forbids the particle to 
travel far from the origin, the size of the localization region being ^/h/ {muj). This so- 
called strong localization has been thoroughly investigated in physics and mathematics 
[31 [3371 [3381 137^14211 [44^ 

7.2. Anderson localization. The previous example of a single quantum har- 
monic well is too idealized. A piece of matter contains an extremely large number of 
interacting atoms. Even if one focuses onto a single atom in an effective potential, 
the form of this potential may be so complicated that the study of the underlying 
eigenfunctions would in general be intractable. In 1958, Anderson considered a lat- 
tice model for a charge carrier in a random potential and proved the localization of 
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(a) (b) ^ 



Fig. 7.1. Illustration of the Anderson transition in a tight-binding model (or so-called SU(2) 
model) in the two-dimensional symplectic class ^20, 21, 372, 373L Three shown eigenf unctions (with 
the energy close to 1) were computed for three disorder strengths W that correspond to (a) metallic 
state (W < Wo), (b) critical state (W = Wq), and (c) insulating state (W > Wq), Wq = 5.952 being 
the critical disorder strength. The latter eigenfunction is strongly localized that prohibits diffusion 
of charge carriers (i.e., no electric current). The eigenf unctions were computed and provided by H. 
Obuse (unpublished earlier). 



eigenfunctions under certain conditions [9 . The localization of charge carriers means 
no electric current through the medium (insulating state), in contrast to metallic or 
conducting state when the charge carriers are not localized. The Anderson transition 
between insulating and conducting states is illustrated for the tight-binding model 



on Fig. 7.1 The shown eigenfunctions were computed by Obuse for three disor- 
der strengths W that correspond to metallic {W < Wq)^ critical {W = Wq)^ and 
insulating (W > Wq) states, I^o = 5.952 being the critical disorder strength. The lat- 
ter eigenfunction is strongly localized that prohibits diffusion of charge carriers (i.e., 
no electric current). The Anderson localization which explains the metal-insulator 
transitions in semiconductors, was thoroughly investigated during the last fifty years 
(see [561 EH ISH [Ml [4281 ME [SSI [lS2] for details and references). Similar 
localization phenomena were observed for microwaves with two-dimensional random 
scattering [142^ , for light in a disordered medium [504^ and in disordered photonic crys- 
tals |44H 1452] , for matter waves in a controlled disorder [70] and in non-interacting 
Bose-Einstein condensate [429 , and for ultrasound [248\ The multifractal structure 
of the eigenfunctions at the critical point (illustrated by Fig. |7.1| 3) has also been 
intensively investigated (see [ 17Q[ 1220] and references therein). Localization of eigen- 
states and transport phenomena in one-dimensional disordered systems are reviewed 
in [251j . An introduction to wave scattering and the related localization is given in 



7.3. Trapping in infinite waveguides. In both previous cases, localization 
of eigenfunctions was related to an external potential. In particular, if the potential 
was not strong enough, Anderson localization could disappear (Fig. |7.1^ ). Is the 
presence of a potential necessary for localization? The formal answer is positive 
because the eigenstates of the Laplace operator in the whole space are simply 
^iik-x) (parameterized by the vector k) which are all extended in R^. These waves are 
called "resonances" (not eigenfunctions) of the Laplace operator, as their I/2-norm is 
infinite. 

The situation is different for the Laplace operator in a bounded domain with 
Dirichlet boundary condition. In quantum mechanics, such a boundary presents a 
"hard wall" that separates the interior of the domain with zero potential from the ex- 
terior of the domain with infinite potential. For instance, this "model" was employed 
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by Crommie et al. to describe the confinement of electrons to quantum corrals on 
a metallic surface |14Q| (see also their figure 2 that shows the experimental spatial 
structure of the electron's wavefunction) . Although the physical interpretation of a 
boundary through an infinite potential is instructive, we will use the mathematical 
terminology and speak about the eigenvalue problem for the Laplace operator in a 
bounded domain without potential. 

For unbounded domains, the spectrum of the Laplace operator consists of two 
parts: (i) the discrete (or point-like) spectrum, with eigenfunctions of finite L2 norm 
that are necessarily "trapped" or "localized" in a bounded region of the waveguide, 
and (ii) the continuous spectrum, with associated functions of infinite L2 norm that 
are extended over the whole domain. The continuous spectrum may also contain 
embedded eigenvalues whose eigenfunctions have finite L2 norm. A wave excited at 
the frequency of the trapped eigenmode remains in the localization region and does 
not propagate. In this case, the definition (7.1) of localization is again rigorous, as 
for any bounded subset VLq of an unbounded domain Vt^ one has jUdi^o) / fJ^di^) = 0, 
while the ratio of L2 norms can be made arbitrarily small by expanding Qq. 

This kind of localization in classical and quantum waveguides has been thoroughly 
investigated (see reviews |163H328j and also references in [376]). In the seminal paper, 
Rellich proved the existence of a localized eigenfunction in a deformed infinite cylinder 
[422 . His results were significantly extended by Jones j264j . Ursell reported on the 
existence of trapped modes in surface water waves in channels [48714489] , while Parker 
observed experimentally the trapped modes in locally perturbed acoustic waveguides 
[385, 386 . Exner and Seba considered an infinite bent strip of smooth curvature 
and showed the existence of trapped modes by reducing the problem to Schrodinger 
operator in the straight strip, with the potential depending on the curvature [171 . 
Goldstone and Jaffe gave the variational proof that the wave equation subject to 
Dirichlet boundary condition always has a localized eigenmode in an infinite tube of 
constant cross-section in any dimension, provided that the tube is not exactly straight 
^06 . This result was further extended by Chenaud et al to arbitrary dimension 
|127j . The problem of localization in acoustic waveguides with Neumann boundary 
condition has also been investigated [167', '168]. For instance, Evans et al. considered a 
straight strip with an inclusion of arbitrary (but symmetric) shape [168 (see [150] for 
further extensions). Such an inclusion obstructed the propagation of waves and was 
shown to result in trapped modes. The effect of mixed Dirichlet, Neumann and Robin 
boundary conditions on the localization was also investigated (see [971 11571 119H 1376] 
and references therein). A mathematical analysis of guided water waves was developed 
by Bonnet-Ben Dhia and Joly [82 (see also \ST). Lower bounds for the eigenvalues 
below the cut-off frequency (for which the associated eigenfunctions are localized) were 
obtained by Ashbaugh and Exner for infinite thin tubes in two and three dimensions 
[22] . In addition, these authors derived an upper bound for the number of the trapped 
modes. Exner et al. considered the Laplacian in finite-length curved tubes of arbitrary 
cross-section, subject to Dirichlet boundary conditions on the cylindrical surface and 
Neumann conditions at the ends of the tube. They expressed a lower bound for 
the spectral threshold of the Laplacian through the lowest eigenvalue of the Dirichlet 
Laplacian in a torus determined by the geometry of the tube |173j . In a different work, 
Exner and co-worker investigated bound states and scattering in quantum waveguides 
coupled laterally through a boundary window jl72j . 

Examples of waveguides with numerous localized states were reported in the lit- 
erature. For instance, Avishai et al. demonstrated the existence of many localized 
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Fig. 7.2. Two examples of a bounded domain 17 = r2i U with a branch of variable cross- 
sectional profile. When the eigenvalue X is smaller than the cut-off 'frequency" ji, the associated 
eigenfunction exponentially decays in the branch Q2 ctnd is thus mainly localized in fli . Note that 
the branch itself may even be increasing. 




Fig. 7.3. Examples of localized Dirichlet Laplacian eigenf unctions with an exponential decay: 
square with a branch (from jl51f }, L-shape and crossing of two stripes (from \152^ ). 

states for a sharp "broken strip", i.e., a waveguide made of two channels of equal 
width intersecting at a small angle [32]. Carini and co-workers reported an ex- 
perimental confirmation of this prediction and its further extensions [116 | (117| I331j . 
Bulgakov et al. considered two straight strips of the same width which cross at an 
angle 9 G (0,7r/2) and showed that, for small 6>, the number of localized states is 
greater than (1 — 2~^/^)^/^/^ [96 . Even for the simple case of two strips crossed at 
the right angle ^ = 7r/2, Schult et al. showed the existence of two localized states, one 
lying below the cut-off frequency and the other being embedded into the continuous 
spectrum |451j . 

7.4. Exponential estimate for eigenfunctions. Qualitatively, an eigenmode 
is trapped when it cannot "squeeze" outside the localization region through narrow 
channels or branches of the waveguide. This happens when typical spatial variations 
of the eigenmode, which are in the order of a wavelength ttA"^/^, are larger than the 
size a of the narrow part, i.e., ttA"^/^ > a or \ < tt'^ / o? [253]. This simplistic argu- 
ment suggests that there exists a threshold value fi (which may eventually be 0), or 
so-called cut-off frequency, such that the eigenmodes with X < ja are localized. More- 
over, this qualitative geometrical interpretation is well adapted for both unbounded 
and bounded domains. While the former case of infinite waveguides was thoroughly 
investigated, the existence of trapped or localized eigenmodes in bounded domains 
has attracted less attention. Even the definition of localization in bounded domains 
remains conventional because all eigenfunctions have finite L2 norm. 

This problem was studied by Delitsyn and co-workers for domains with branches 
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n = 1 



n = 2 



n = 3 



Fig. 7.4. T/ie /irst t/iree Dirichlet Laplacian eigenf unctions for three elongated domains: (a) 
rectangle of size 25 x 1; (b) right trapezoid with bases 1 and 0.9 and height 25 which is very close to 
the above rectangle; and (c) right triangle with edges 25 and 1 (half of the rectangle). There is no 
localization for the first shape, while the first eigenf unctions for the second and third domains tend 
to be localized. 

of variable cross-sectional profiles [151 . More precisely, one considers a bounded 
domain Q C R"^ {d = 2,3,...) with a piecewise smooth boundary dft and denote 
Q{z) = fl n {x G R^ : xi = z} the cross-section of at Xi = 2: G M by a hyperplane 
perpendicular to the coordinate axis xi (Fig. 7.2). Let 



zi = mf{z e R : Q{z) ^ 0}, Z2 = sup{z G R : Q{z) ^ 0}, 

and we fix some zq such that Zi < zq < Z2. Let be the first eigenvalue of 

the Laplace operator in Q{z), with Dirichlet boundary condition on dQ{z), and ja = 
inf /i(^). Let be a Dirichlet-Laplacian eigenfunction in and A the associate 

eigenvalue. If A < /i, then 



\ML2iQiz)) < II^IIl2(Q(^o)) exp(-/3v'/i- a {z - zo)) {z > zo), (7.6) 

with (3 = 1/a/2. Moreover, if (ei • n{x)) > for all x G dft with xi > zq, where ei is 
the unit vector (1, 0, 0) in the direction xi, and n{x) is the normal vector at x G dCl 
directed outwards the domain, then the above inequality holds with {3 = 1. 

In this statement, a domain ^ is arbitrarily split into two subdomains, l^i (with 
xi < zq) and 0.2 (with xi > zq), by the hyperplane at xi = zq (the coordinate axis xi 
can be replaced by any straight line). Under the condition \ < ji^ the eigenfunction 
u exponentially decays in the subdomain 0.2 which is loosely called "branch". Note 
that the choice of the splitting hyperplane (i.e., 2:0) determines the threshold ji. 
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The theorem formahzes the notion of the cut-off frequency /i for branches of 
variable cross-sectional profiles and provides a constructive way for its computation. 
For instance, if Q2 is a rectangular channel of width a, the first eigenvalue in all cross- 
sections Q{z) is TT^/a^ (independently of z) so that fi = 7r^/a^, as expected. The 
exponential estimate quantifies the "difficulty" of penetration, or "squeezing", into 
the branch Q2 and ensures the localization of the eigenfunction u in Qi. Since the 
cut-off frequency /i is independent of the subdomain r^i, one can impose any boundary 
condition on dCli (that still ensures the self- adjoint ness of the Laplace operator). In 
turn, the Dirichlet boundary condition on the boundary of the branch Q2 is relevant, 
although some extensions were discussed in |151j . It is worth noting that the theorem 
also applies to infinite branches 1^2, under supplementary condition /i(z) ^ 00 to 
ensure the existence of the discrete spectrum. 

According to this theorem, the L2-norm of an eigenfunction with X < jii in Q{z) = 
Qr]{x : xi > z} can be made exponentially small provided that the branch Q2 



is long enough. Taking Qq = Q\Q{z), the ratio of L2-norms in Eq. (7.1) can be made 
arbitrarily small. However, the second ratio may not be necessarily small. In fact, its 
smallness depends on the shape of the domain Q. This is once again a manifestation 
of the conventional character of localization in bounded domains. 



Figure 7.3 presents several examples of localized Dirichlet Laplacian eigenfunc- 
tions showing an exponential decay along the branches. Since an increase of branches 
diminishes the eigenvalue A and thus further enhances the localization, the area of the 
localized region Qi can be made arbitrarily small with respect to the total area (one 
can even consider infinite branches). Examples of an L-shape and a cross illustrate 
that the linear sizes of the localized region do not need to be large in comparison 
with the branch width (a sufficient condition for getting this kind of localization was 
proposed in |152j ). It is worth noting that the separation into the localized region and 



branches is arbitrary. For instance. Fig. 7.4 shows several localized eigenfunctions for 



elongated triangle and trapezoid, for which there is no explicit separation. 

Localization and exponential decay of Laplacian eigenfunctions were observed for 
various perturbations of cylindrical domains [ 115[ 1272] 13651 . For instance, Kamotskii 
and Nazarov studied localization of eigenfunctions in a neighborhood of the lateral 
surface of a thin domain |272j . Nazarov and co-workers analyzed the behavior of 



eigenfunctions for thin cylinders with distorted ends |115H365j . For a bounded domain 
uj C R^~^ (n > 2) with a simple closed Lipschitz contour duo and Lipschitz functions 
H±{r]) in cc; = ujUduj^ the thin cylinder with distorted ends is defined for a given small 
£ > as 

= {{x,y) eRxW-^ : y = er], - 1 - eH_{r]) < x < 1 ^ eH^{r]), r] e uj}. 
One can view this domain as a thin cylinder [—1,1] x (eu) to which two distorted 



"cups" characterized by functions H±^ are attached (Fig. 7.5 1). The Neumann bound- 
ary condition is imposed on the curved ends F^: 

F^ = {(x, y) eRx : y = er], x = ±1 ± eH±{r]), r] e u}, 

while the Dirichlet boundary condition is set on the remaining lateral side of the 
domain: = dQ^\{r^_^ U F^). When the ends of the cylinder are straight {H± = 0), 
the eigenfunctions are factored as Umn{x,y) = cos(7rm(x + l)/2)(pn{y)^ where (pn{y) 
are the eigenfunctions of the Laplace operator in the cross-section uj with Dirichlet 
boundary condition. These eigenfunctions are extended over the whole cylinder, due 
to the cosine factor. Nazarov and co-workers showed that distortion of the ends (i.e.. 
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(a) 



Dirichlet 



Dirichlet 



(b) 



-H . 



Dirichlet 



Dirichlet 



(c) (d) (e) (f) 




Fig. 7.5. (a) A thin cylinder fl^ G of width e with two distorted ends defined by functions 
H±{rf), ?7 G a; = (0, 1). (b) In the limit e ^ 0, the analysis is reduced to a semi-infinite cylinder VL 
with one distorted end. (c,d,e,f) Four semi-infinite cylinders with various distorted ends (top), the 
first eigenfunction for the Laplace operator in these shapes with mixed Dirichlet- Neumann boundary 
condition (middle), and the first eigenfunction for the Laplace o perator with purely Dirichlet bound- 
ary condition. According to the sufficient conditions \7.'1\7.8\l , the first eigenfunction is localized 
near the distorted end for cases ^d\ ^e^ and f\ and not localized for the case ^c\ No localization 
happens when the Dirichlet boundary condition is set over the whole domain. For numerical compu- 
tation, semi-infinite cylinders were Hruncated" and auxilary Dirichlet boundary condition was set 
at the right straight end. 

H± ^ 0) may lead to localization of the ground eigenfunction in one (or both) ends, 
with an exponential decay toward the central part. In the limit £ ^ 0, the thinning 
of the cylinder can be alternatively seen as its outstretching, allowing one to reduce 
the analysis to a semi-infinite cylinder with one distorted end (Fig. |7.5| 3), described 
by a single function H{r]): 

n = {(x, T]) eRx : - H{r]) < x, r] e u}. 

Two sufficient conditions for getting the localized ground eigenfunction at the dis- 
torted end were proposed in [115 : 

(i) For H G C(cj), the sufficient condition reads 

y dT^ H{7^)(^\vMv)\' - <o, (7.7) 

where /ii is the smallest eigenvalue corresponding to cpi in the cross-section uj (in two 
dimensions, when uj = [—1/2,1/2], one has (pi{r]) = sin7r(7^ + 1/2) and /ii = tt^ so 
that this condition reads Xl(/2 dv H(r]) cos(27r7^) > |365j ): 
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(ii) Under a stronger assumption H G the sufficient condition simplifies to 

j dr,[^i{il)f^H{r,)<0. (7.8) 

This inequality becomes true for a subharmonic profile —H (i.e., for AH{rj) < 0) with 
but false for super harmonicj^ 

Figure [73] shows several examples for which the sufficient condition is either sat- 



isfied (7.5i,e,f), or not ( |7.5| 3). Nazarov and co-workers showed that these results are 
applicable to bounded thin cylinders for small enough e. In addition, they found 
out a domain where the first eigenfunction concentrates at the both ends simulta- 
neously. Finally, they showed that no localization happens in the case in which the 
mixed Dirichlet-Neumann boundary condition is replaced by the Dirichlet boundary 



condition onto the whole boundary, as illustrated on Fig. 7.5 (see |115j for further 
discussions and results). 

Friedlander and Solomyak studied the spectrum of the Dirichlet Laplacian in a 
family of narrow strips of variable profile: = {(x,y) G : — a < x < 6, < 
y < £h{x)} \195\ 1196] . The main assumption was that x = is the only point of 
global maximum of the positive, continuous function h{x). In the limit £ ^ 0, they 
found the two-term asymptotics of the eigenvalues and the one-term asymptotics 
of the corresponding eigenfunctions. The asymptotic formulas obtained involve the 
eigenvalues and eigenfunctions of an auxiliary ODE on R that depends only on the 
behavior of h{x) as x ^ 0, i.e., in the vicinity of the widest cross-section of the strip. 

7.5. Dumbbell domains. Yet another type of localization emerges for domains 
that can be split into two or several subdomains with narrow connections (of "width" 
s) |417] . a standard example being a dumbbell: Cl^ = Qi U U Q2 (Fig- |7.6^ ). 
The asymptotic behavior of eigenvalues and eigenfunctions in the limit e ^ was 
thoroughly investigated for both Dirichlet and Neumann boundary conditions [261 . 
We start by considering Dirichlet boundary condition. 

In the limiting case of zero width connections, the subdomains Qi (z = 1, A/") 
become disconnected, and the eigenvalue problem can be independently formulated 
for each subdomain. Let be the set of eigenvalues for the subdomain Qi. Each 
Dirichlet eigenvalue of the Laplace operator in the domain approaches to an 
eigenvalue corresponding to one limiting subdomain Qi C ft^: G for certain 
i. Moreover, if 



A,nA, =0 Vi^j, (7.9) 

the space of eigenfunctions in the limiting (disconnected) domain is the direct 
product of spaces of eigenfunctions for each subdomain fti (see [143 for discussion 
on convergence and related issues). This is a basis for what we will call "bottle-neck 
localization". In fact, each eigenfunction on the domain approaches an eigen- 
function of the limiting domain which is fully localized in one subdomain Qi and 
zero in the others. For a small s, the eigenfunction is therefore mainly localized in 
the corresponding i-th subdomain 1]^, and is almost zero in the other subdomains. In 
other words, for any eigenfunction, one can take the width e small enough to ensure 
that the L2-norm of the eigenfunction in the subdomain Qi is arbitrarily close to that 



^ In the discussion after Theorem 3 in [115) . the sign minus in front of H was omitted. 
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Fig. 7.6. (a) A dumbbell domain fl^ is the union of two bounded domains fli, fl2 ctnd a 
narrow ^'connector'' of width e. (b) In the limit e ^ 0, the connector degenerates to a curve 
(here, an interval) so that the subdomains Vli and 0.2 become disconnected, (c) In Beale's work, 
is composed of two components, a bounded domain Qi and unbounded domain Q2, which are 
connected by a narrow channel . 



n = 1 n = 7 n = 8 n = 11 




Fig. 7.7. Several Dirichlet Laplacian eigenf unctions for a dumbbell domain which is composed 
of two rectangles and connected by the third rectangle (from jl51f }. The 1st and 7th eigenf unctions 
are localized in the larger subdomain, the 8th eigenfunction is localized in the smaller subdomain, 
while the 11th eigenfunction is not localized at all. Note that the width of connection is not small 
(1/A of the width of both subdomains) . 



in the whole domain 



Vm>l 3i G {1,...,7V} V(5g(0,1) 3£>0 



>(1 



(7.10) 

This behavior is exemphfied for a dumbbeh domain which is composed of two rectan- 
gles and connected by the third rectangle (Fig. 7.7). The 1st and 7th eigenfunctions 
are localized in the larger rectangle, the 8th eigenfunction is localized in the smaller 
rectangle, while the 11th eigenfunction is not localized at all. Note that the width of 
connection is not too small (1/4 of the width of both subdomains). 

It is worth noting that, for a small fixed width e and a small fixed threshold ^, 
there may be infinitely many high-frequency "non-localized" eigenfunctions, for which 
the inequality ( 7.10 ) is not satisfied. In other words, for a given connected domain with 
a narrow connection, one can only expect to observe a finite number of low-frequency 
localized eigenfunctions. We note that the condition (7.9) is important to ensure that 
limiting eigenfunctions are fully localized in their respective subdomains. Without 
this condition, a limiting eigenfunction may be a linear combination of eigenfunctions 
in different subdomains with the same eigenvalue that would destroy localization. 
Note that the asymptotic behavior of eigenfunctions at the "junction" was studied by 
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Felli and Terracini |178j . 

For Neumann boundary condition, the situation is more complicated, as the eigen- 
values and eigenfunctions may also approach the eigenvalues and eigenfunctions of the 
limiting connector (in the simplest case, the interval). Arrieta considered a planar 
dumbbell domain consisted of two disjoint domains fti and connected by a 
channel of variable profile g{x): = {x eM? : < xi < 1, < X2 < 
where g G C^(0, 1) and g{xi) > for all xi G [0, 1]. In the limit £ ^ 0, each eigen- 
value of the Laplace operator in with Neumann boundary condition was shown to 
converge either to an eigenvalue jj^k of the Neumann-Laplace operator in Qi U 1^2, or 
to an eigenvalue of the Sturm-Liouville operator —^{gux)x acting on a function 
u on (0,1), with Dirichlet boundary condition [l4l|T5]. The first-order term in the 
small e-asymptotic expansion was obtained. The special case of cylindrical channels 
(of constant profile) in higher dimensions was studied by Jimbo [258 (see also results 
by Hempel et al. |237j ). Jimbo and Morita studied an A^-dumbbell domain, i.e., a 
family of N pairwise disjoint domains connected by thin channels [259 . They proved 
that = Cm£^~^ + o{e^~^) as £ ^ for m = 1, 2, . . . , A^, while A^^^^^ uniformly 
bounded away from zero, where d is the dimension of the embedding space, and Cm 
are shape-dependent constants. Jimbo also analyzed the asymptotic behavior of the 
eigenvalues A^ with m > N under the condition that the sets {fik} and {vk} do not 
intersect |26Qj . In particular, for an eigenvalue A^ that converges to an element of 
{/i/c}, the asymptotic behavior is A^ = /i^ + Cm£^~^ + o{e^~^). 

Brown and co-workers studied upper bounds for |A^ — A^| and showed ^89j: 

(i) If Aie{Mfc}\K}, 

|A^-A^|<C|ln£|-i/2 (d=2), 

(ii) If A^ e {uk} \ {^ik}, 

|AL-A^|<Cei/2|ln£| {d = 2), 
|A^-A^|<Csi/2 (d>3). 

For a dumbbell domain in with a thin cylindrical channel of a smooth profile, 
Gadyl'shin obtained the complete small e asymptotics of the Neumann-Laplace eigen- 
values and eigenfunctions and explicit formulas for the first term of these asymptotics, 
including multiplicities [198r.2QQj . 

Arrieta and Krejcirik considered the problem of spectral convergence from an- 
other point of view [19]. They showed that if C fls are bounded domains and 
if the eigenvalues and eigenfunctions of the Laplace operator with Neumann bound- 
ary condition in converge to the ones in then necessarily fid{^s\^o) ^ as 
£ ^ 0, while it is not necessarily true that dist(l]£, Qq) 0. As a matter of fact, they 
constructed an example of a perturbation where the spectra behave continuously but 
dist(r^£, Qq) ^ 00 as £ ^ 0. 

A somewhat related problem of scattering frequencies of the wave equation as- 
sociated to an exterior domain in with an appropriate boundary condition was 
investigated by Beale [55] (for more general aspects of geometric scattering theory, 
see [355 ). We recall that a scattering frequency >/A of an unbounded domain Q is 
a (complex) number for which there exists a nontrivial solution of Au -\- Xu = in 
subject to Dirichlet, Neumann, or Robin boundary condition and to an "outgo- 
ing" condition at infinity. In Beale's work, a bounded cavity Qi was connected by a 
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thin channel to the exterior (unbounded) space 1^2- More specificahy, he considered a 
bounded domain D such that its complement in has a bounded component Vti and 
an unbounded component 0.2- After that, a thin "hole" in D was made to connect 
both components (Fig. |7.6| 3). Beale showed that the joint domain = U U 1^2 
with Dirichlet boundary condition has a scattering frequency which is arbitrarily close 
either to an eigenfrequency (i.e., the square root of the eigenvalue) of the Laplace op- 
erator in r^i, or to a scattering frequency in 0.2^ provided the channel is narrow 
enough. The same result was extended to Robin boundary condition of the form 
du / dn -\- hu = on dCt^ ^ where is a function on dCt^ with a positive lower bound. 
In both cases, the method in his proof relies on the fact that the lowest eigenvalue 
of the channel tends to infinity as the channel narrows. However, it is no longer 
true for Neumann boundary condition. In this case, with some restrictions on the 
shape of the channel, Beale proved that the scattering frequencies converge not only 
to the eigenfrequencies of Qi and scattering frequencies of Q2 but also to the longitu- 
dinal frequencies of the channel. Similar results can be obtained in domains of space 
dimension other than 3. 

There are other problems for partial differential equations in dumbbell domains 
which undergo a singular perturbation [225' "240' I341j . For instance, in a series of 
articles, Arrieta et al. studied the behavior of the asymptotic nonlinear dynamics 
of a reaction-diffusion equation with Neumann boundary condition p^HT8| . In this 
context, dumbbell domains appear naturally as the counterpart of convex domains for 
which the stable stationary solutions to a reaction-diffusion equation are necessarily 
spatially constant |123j . As explained in [16 , one way to produce "patterns" (i.e., 
stable stationary solutions which are not spatially constant), is to consider domains 
which make it difficult to diffuse from one part of the domain to the other, making 
a constriction in the domain. Kosugi studied the convergence of the solution to a 
semilinear elliptic equation in a thin network-shaped domain which degenerates into 
a geometric graph when a certain parameter tends to zero [288 (see also [ 295^ I4Q7[ 

WW)- 

7.6. Localization in irregularly-shaped domains. As we have seen, a nar- 
row connection between subdomains could lead to localization. How narrow should 
it be? A rigorous answer to this question is only known for several "tractable" cases 



such as dumbbell-like or cylindrical domains (Sec. 7.5). Sapoval and co-workers have 
formulated and studied the problem of localization in irregularly-shaped or fractal 
domains through numerical simulations and experiments [1691 I175[ I224[ I232[ I433[ 
I434[ I442H444] . In the first publication, they monitored the vibrations of a prefractal 
"drum" (i.e., a thin membrane with a fixed boundary) which was excited at different 
frequencies [442]. Tuning the frequency allowed them to directly visualize different 
Dirichlet Laplacian eigenfunctions in a (prefractal) quadratic von Koch snowflake (an 



example is shown on Fig. 7.8). For this and similar domains, certain eigenfunctions 



were found to be localized in a small region of the domain, for both Dirichlet and 



Neumann boundary conditions (Fig. 7.9). This effect was first attributed to self- 



similar structure of the domain. However, similar effects were later observed through 



numerical simulations for non-fractal domains [175, 445^ . as illustrated by Fig. 7.10 
In the study of sound attenuation by noise-protective walls, Felix and co-workers have 
further extended the analysis to the union of two domains with different refraction 
indices which are separated by an irregular boundary [175. 176. 445^ . Many eigen- 
functions of the related second order elliptic operator were shown to be localized on 
this boundary (so-called "astride localization"). A rigorous mathematical theory of 
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^7 = l 




Fig. 7.8. The unit square and three prefractal domains obtained iteratively one from the 
other (two sides of these domains are finite generations of the Von Koch curve of fractal dimension 
S/2). These domains were intensively studied, both numerically and experimentally, by Sapoval and 
co-workers fl69l \224\ \232\ \433\ \434\ \442W444j . 
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Fig . 7.9 . Several Dirichlet (top) and Neumann (bottom) eigenf unctions for the third domain 
on Fig. \7.8\ (g = 2). The 38th Dirichlet and the 12th N euma nn eigenf unctions are localized in a 
small subdomain (located in the upper right corner on Fig. \ 7. 8}^ , while the first/second Dirichlet and 
the Ath Neumann eigenf unctions are almost zero on this subdomain. Finally, the 8th Dirichlet and 
the second Neumann eigenf unctions are examples of eigenf unctions extended over the whole domain. 



these important phenomena is stih missing. Takeda al. observed experimentahy the 
electromagnetic field at specific frequency to be confined in the central part of the 
third stage of three-dimensional fractals called the Menger sponge [480 . This local- 
ization was attributed to a singular photon density of states realized in the fractal 
structure. 

A number of mathematical studies were devoted to the theory of partial differen- 
tial equations on fractals in general and to localization of Laplacian eigenfunctions in 
particular (see [28Ql |477] and references therein). For instance, the spectral proper- 
ties of the Laplace operator on Sierpinski gasket and its extensions were thoroughly 
investigated j42H44| \46\ \75\ \197\ I456j . Barlow and Kigami studied the localized eigen- 
functions of the Laplacian on a more general class of self-similar sets (so-called post 
critically finite self-similar sets, see [281, 282 for details). They related the asymptotic 
behavior of the eigenvalue counting function to the existence of localized eigenfunc- 
tions and established a number of sufficient conditions for the existence of a localized 
eigenfunction in terms of the symmetries of a set [45 . 

Berry and co-workers developed a new method to approximate the Neumann 
spectrum of the Laplacian on a planar fractal set ^7 as a renormalized limit of the 
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Fig. 7.10. Examples of localized Neumann Laplacian eigenf unctions in two domains adapted 
from I175L- square with many elongated holes (top) and random sawteeth (bottom). Colors represent 
the amplitude of eigenf unctions, from the most negative value (dark blue), through zero (green), 
to the largest positive value (dark red). One can notice that the eigenf unctions on the top are 
not negligible outside the localization region. This is yet another illustration for the conventional 
character of localization in bounded domains. 



Neumann spectra of the standard Laplacian on a sequence of domains that approxi- 
mate Q from the outside |66J. They apphed this method to compute the Neumann- 
Laplacian eigenfunctions in several domains, including a sawtooth domain, Sierpinski 
gasket and carpet, as well as nonsymmetric and random carpets and the octagasket. 
In particular, they gave a numerical evidence for the localized eigenfunctions for a 
sawtooth domain, in agreement with the earlier work by Felix et al. [175^ . 

Heilman and Strichartz reported several numerical examples of localized Neumann- 



Laplacian eigenfunctions in two domains [233 , one of them is illustrated on Fig. 7.11 i. 
Each of these domains consists of two subdomains with a narrow but not too narrow 
connection. This is a kind of dumbbell shape with a connector of zero length. Heil- 
man and Strichartz argued that one subdomain must possess an axis of symmetry 
for getting localized eigenfunctions. Since an anti-symmetric eigenfunction vanishes 
on the axis of symmetry, it is necessarily small near the bottle-neck that somehow 
"prevents" its extension to the other domain. Although the argument is plausible, 
we have to stress that such a symmetry is neither sufficient, nor necessary for lo- 
calization. It is obviously not sufficient because even for symmetric domain, there 
exist plenty of extended eigenfunctions (including the trivial example of the ground 
eigenmode which is a constant over the whole domain). In order to illustrate that the 
reflection symmetry is not necessary, we plot on Fig. |7.11| 3,c examples of localized 
eigenfunctions for modified domains for which the symmetry is broken. Although 
rendering the upper domain less and less symmetric gradually reduces or even fully 
destroys localization (Fig. 7.111), its "mechanism" remains poorly understood. We 



also note that methods of Sec. |7.4| are not applicable in this case because of Neumann 
boundary condition. 

Lapidus and Pang studied the boundary behavior of the Dirichlet Laplacian eigen- 
functions and their gradients on a class of planar domains with fractal boundary, 
including the triangular and square von Koch snowflakes and their polygonal approx- 
imations |3Q5j . A numerical evidence for the boundary behavior of eigenfunctions 



was reported in [307 , with numerous pictures of eigenfunctions. Later, Daubert and 
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Fig. 7.11. Neumann- Laplace eigenfunction U4 in the original "cow" domain from f233f (a) 
and in three modified domains (h,c,d), in which the reflection symmetry of the upper subdomain 
is broken. The fourth eigenfunction is localized for the first three domains (a,b,c), while the last 
domain with the strongest modification shows no localization (d). Colors represent the amplitude 
of an eigenfunction, from the most negative value (dark blue), through zero (green), to the largest 
positive value (dark red). 



Lapidus considered more specifically the localization character of eigenfunctions in 
von Koch domains [145 . In particular, different "measures" of localization were dis- 
cussed. 

Note also that Filoche and Mayboroda studied the problem of localization for 
bi-Laplacian in rigid thin plates and discovered that clamping just one point inside 
such a plate not only perturbs its spectral properties, but essentially divides the plate 
into two independently vibrating regions |185j . 

7.7. High-frequency localization. A hundred years ago. Lord Rayleigh doc- 
umented an interesting acoustical phenomenon in the whispering gallery under the 
dome of Saint Paul's Cathedral in London ^419^ (see also [415" 416 ). A whisper 
of one person propagated along the curved wall to another person stood near the 
wall. Keller and Rubinow discussed the related "whispering gallery modes" and also 
"bouncing ball modes", and showed that these modes exist for a two-dimensional 
domain with arbitrary smooth convex curve as its boundary j277j . A semiclassi- 
cal approximation of Laplacian eigenfunctions in convex domains was developed by 
Lazutkin [33l[3TT^3TT (see also [11 [4131 SH [467]). Chen and co-workers analyzed 
Mathieu and modified Mathieu functions and reported another type of localization 
named "focusing modes" [126 . All these eigenmodes become more and more localized 
in a small subdomain when the associated eigenvalue increases. This so-called high- 
frequency or high-energy limit was intensively studied for various domains, named as 
quantum billiards |222[ I235| I254[ [447| 1474] . In quantum mechanics, this limit is known 
as semi-classical approximation [61]. In optics, it corresponds to ray approximation 
of wave propagation, from which the properties of an optical, acoustical or quantum 
system can often be reduced to the study of rays in classical billiards. Jakobson et 
al. gave an overview of many results on geometric properties of the Laplacian eigen- 
functions on Riemannian manifolds, with a special emphasis on high-frequency limit 
(weak star limits, the rate of growth of Lp norms, relationships between positive and 
negative parts of eigenfunctions, etc.) [25T (see also [4*, "486]). Bearing in mind this 
comprehensive review, we start by illustrating the high-frequency localization and 
the related problems in simple domains such as disk, ellipse and rectangle for which 
explicit estimates are available. After that, some results for quantum billiards are 
summarized. 



7.7.1. Whispering gallery and focusing modes. The disk is the simplest 
shape for illustrating the whispering gallery and focusing modes. The explicit form 
(3.9) of eigenfunctions allows one to get accurate estimates and bounds, as shown 
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below. When the index k is fixed and n increases, the Bessel functions Jn{oink'f / R) 
become strongly attenuated near the origin (as Jn(^) ~ {z/2)'^ /n\ at small z) and 
essentially localized near the boundary, yielding whispering gallery modes. In turn, 
when n is fixed and k increases, the Bessel functions rapidly oscillate, the amplitude 
of oscillations decreasing towards the boundary. In that case, the eigenfunctions are 
mainly localized at the origin, yielding focusing modes. 

These qualitative arguments were rigorously formulated in j37Qj . For each eigen- 
function Unk on the unit disk 1^, one introduces the subdomain Vtrik = {x ^M? : \x\ < 
dn/(^nk} C where dn = n — n^/^, and ank are, depending on boundary conditions, 
the positive zeros of either Jn{z) (Dirichlet), or J^iz) (Neumann) or J^i^) + hJn{z) 
for some > (Robin), with n = 0, 1, 2, ... denoting the order of Bessel function Jn{z) 
and k = 1, 2, 3, ... counting zeros. Then for any p>l (including p = oo), there exists 
a universal constant Cp > such that for any k = 1, 2, 3, ... and any large enough n, 
the Laplacian eigenfunction Unk for Dirichlet, Neumann or Robin boundary condition 
satisfies 

"^^"^^^-^^ < c,nH* exp(-nV3 ln(2)/3). (7.11) 

WunkUpin) 

The definition of Qnk and the above estimate imply 

hm f = 0, hm --— = 1. (7.12) 

This theorem shows the existence of infinitely many Laplacian eigenmodes which are 
Lp-localized in a thin layer near the boundary dVt. In fact, for any prescribed thresh- 



olds for both ratios in (7.1 ), there exists no such that for all n > no, the eigenfunctions 



Ink are Lp-localized. These eigenfunctions are called "whispering gallery eigenmodes" 



and illustrated on Fig. 7.12 



We outline a peculiar relation between high-frequency and low-frequency localiza- 



tion. The explicit form (3.9) of Dirichlet Laplacian eigenfunctions u^k leads to their 
simple nodal structure which is formed by 2n radial nodal lines and k — 1 circular 
nodal lines. The radial nodal lines split the disk into 2n circular sectors with Dirichlet 
boundary conditions. As a consequence, whispering gallery eigenmodes in the disk 



and the underlying exponential estimate (7.11) turn out to be related to the expo- 



nential decay of eigenfunctions in domains with branches (Sec. 7.4), as illustrated on 
Fig. |7.4| for elongated triangles. 

A simple consequence of the above theorem is that for any p > 1 and any open 
subset V compactly included in the unit disk Vt (i.e., V D dVt = 0), one has 

hm = 0, (7.13) 

n^oo \\Unk\\LpiQ) 

and 



Qualitatively, for any subset there exists a sequence of eigenfunctions that pro- 
gressively "escape" V. 

The localization of focusing modes at the origin is revealed in the limit k ^ oo. 
For each R G (0, 1), one defines an annulus Qr = {x G : R < \x\ < 1} C Q of 
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Fig. 7.12. Formation of whispering gallery modes for the unit disk with Dirichlet boundary 
condition: for a fixed k (k = 1 for top figures and k = 2 for bottom figures), an increase of the index 
n leads to stronger localization of eigenf unctions near the boundary. 



the unit disk Q.. Then, for any n = 0,1,2,..., the Laplacian eigenfunction Unk with 
Dirichlet, Neumann or Robin boundary condition satisfies 



IKkh^jn^) ^ I (1 - (1 < p < 4), 

\\u„k\\L,in) \ (p>4). 



(7.15) 



When the index k increases (with fixed n), the eigenfunctions Unk become more and 
more Lp-locahzed near the origin when p > 4 |37Qj . These eigenfunctions are cahed 
"focusing eigenmodes" and ihustrated on Fig. |7.13| The theorem shows that the 
definition of locahzation is sensitive to the norm: the above focusing modes are Lp- 
locahzed for p > 4, but they are not L^-locahzed for p < 4. Similar results for 
whispering gallery and focusing modes hold for a ball in three dimensions [37Qj • 

7.7.2. Bouncing ball modes. Filled ellipses and elliptical annuli are simple 
domains for illustrating bouncing ball modes. For fixed foci (i.e., a given parameter a 
in the elliptic coordinates in Eq. (3.15 )), these domains are characterized by two radii, 
Rq (i^Q = for filled ellipses) and i?, as O = Ur,0) : Ro<r<R, 0<0< 27r}, while 
the eigenfunctions Unki were defined in Sec. |3.4l For each a G (O, |), we consider an 
elliptical sector inside an elliptical domain il (Fig. |3.1[ ) 

= {{r, 0) : Rq < r < R, e {a,7i - a)l} {n ^ a,2ii - a)} . 

For any p > 1, there exists Ac^^^ > such that for any eigenvalue Xnki > ^a,n^ the 
corresponding eigenfunction Unki satisfies [370 (see also ^68j ) 



ll^n/cdl 



Lp{Q\Qo.) 



\\Unkl\\ 



16a 



a/2 



i/p 



exp (^-ay^Xnki sin Q + 



(7.16) 



where 
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Fig. 7.13. Formation of focusing modes for the unit disk with Dirichlet boundary condition: 
for a fixed n (n = for top figures and n = 1 for bottom figures), an increase of the index k leads 
to stronger localization of eigenf unctions at the origin. 
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Fig. 7.14. Formation of bouncing ball modes Unkl in a filled ellipse of radius R = 1 (top) 
and an elliptical annulus of radii Rq = 0.5 and R = 1 (bottom), with the focal distance a = 1, and 
Dirichlet boundary condition. For fixed n = 1 and I = 1, an increase of the index k leads to stronger 
localization of the eigenfunction near the vertical semi-axis (from \370\ ). 



Given that \nki ^ oo as /c increases (for any fixed n and /), while the area of VLcx can 
be made arbitrarily small by sending a 7r/2, the estimate implies that there are 
infinitely many eigenfunctions Unki which are L^-localized in the elliptical sector Q^- 



k^oo \\Unkl\\Lp{n) 



(7.17) 



These eigenfunctions, which are localized near the minor axis, are called "bouncing 



ball modes" and illustrated on Fig. |7.14[ The above estimate allows us to illustrate 
bouncing ball modes which are known to emerge for any convex planar domain with 
smooth boundary [126| 1277] . At the same time, the estimate is as well applicable to 
elliptical annuli, providing thus an example of bouncing ball modes for non-convex 
domains. 
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7.7.3. Domains without localization. The analysis of geometrical properties 
of eigenfunctions in rectangle-like domains Q = (0,^i) x ... x (O^id) C (with sizes 
^1 > 0, -^c^ > 0) iiiay seem to be the simplest case because the eigenfunctions are 
expressed through sines (Dirichlet) and cosines (Neumann), as discussed in Sec. |3.1[ 



The situation is indeed elementary when all eigenvalues are simple, i.e., (^i/^j)^ are 
not rational numbers for all i ^ j. For any p > 1 and any open subset V C 1^, one 
can prove that [370 

C,{V)^ inf 0, (7.18) 



This property is in sharp contrast to Eq. (7.14) for eigenfunctions in the unit disk (or 



ball). The fact that Cp{V) > for any open subset V means that there is no eigen- 
function that could fully "avoid" any location inside the domain, i.e., there is no Lp- 
localized eigenfunction. Since the set of rational numbers has zero Lebesgue measure, 
there is no L^-localized eigenfunctions in almost any randomly chosen rectangle-like 
domain. 

When at least one ratio {ii/ij)'^ is rational, certain eigenvalues are degenerate, and 
the associated eigenfunctions are linear combinations of products of sines or cosines 



(see Sec. 3.1). Although the computation is still elementary for each eigenfunction. 



it is unknown whether the infimum Cp{V) from Eq. (7.18) is strictly positive or not, 
for arbitrary rectangle-like domain ft and any open subset V. For instance, the most 
general known result for a rectangle ft = (0,^i) x (0,^2) states that C2{V) > for 
any V C Q the form V = (0, ^1) x cj, where uj is any open subset of (0, ^2) |11Q| (see 
also |335j ). Even for the unit square, the statement Cp{V) > for any open subset 
V appears as an open problem. More generally, one may wonder whether Cp{V) is 
strictly positive or not for any open subset V in polygonal (convex) domains. 

7.7.4. Quantum billiards. The above examples of whispering gallery or bounc- 
ing ball modes illustrate that certain high-frequency eigenfunctions tend to be local- 
ized in specific regions of circular and elliptical domains. But what is the structure of 
a high-frequency eigenfunction in a general domain? What are these specific regions 
on which a sequence of eigenfunctions may be localized, and whether do localized 
eigenfunctions exist for a given domain? Answers to these and other related ques- 
tions can be found by relating the high-frequency behavior of a quantum system (in 
our case, the structure of Laplacian eigenfunctions) to the classical dynamics in a 
billiard of the same shape [11 [I29l [2211 [HIl [327]. This relation is also known as 
a semi-classical approximation in quantum mechanics and a ray approximation of 
wave propagation in optics, while the correspondence between classical and quantum 
systems can be shown by the WKB method, Eikonal equation or Feynman path in- 
tegrals |179[ I18Q( 1361] . For instance, the dynamics of a particle in a classical billiard 
is translated into quantum mechanism through the stationary Schrodinger equation 
Hun{x) = EnUn{x) with Dirichlct boundary condition, where the free Hamiltonian 
is H = p^/{2m) = — fi^A/(2m), and the energy is related to the corresponding 
Laplacian eigenvalue = 2mEn/h'^. Since \un{x)\'^ is the probability density for 
finding a quantum particle at x, this density should resemble some classical trajec- 
tory of that particle in the (high-frequency) semi-classical limit {h ^ or m ^ 00). 
In particular, some orbits of a particle moving in a classical billiard may appear 
as "scars" in the spatial structure of eigenfunctions in the related quantum billiard 
[SD M [Sa [2221 [2231 UMl |235l HH IMl im Ea [111 HZH. This effect is 



illustrated on Fig. 7.15| by Liu and co-workers who investigated the localization of 
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Fig. 7.15. Several eigenstates with localization on periodic orbits for the spiral-shaped billiard 
with e = 0.1, from ]330f (Reprinted Fig. 6 with permission from Liu et al, Physical Reviews E, 74, 
046214, (2006). Copyright (2006) by the American Physical Society). 




Fig. 7.16. Examples of chaotic billiards: (a) Bunimovich stadium (union of a square and two 
half - disks ) \M \M UM \^34 , 483, 484I, (b) Sinai's billiard ;463\\f6^, (c) mushroom billiard 
\48\ \100f . and (d) hyperbolic billiard ^2L Many other examples are given in f99f . 



Dirichlet Laplacian eigenfunctions on classical periodic orbits in a spiral-shaped bil- 
liard [330] (see also [3T7] ). 

In the classical dynamics, one may distinguish the domains with regular, inte- 
grable and chaotic dynamics. In particular, for a bounded domain ft with an ergodic 
billiard flow [462], Shnirelman's theorem (also known as quantum ergodicity theorem 
jl36[ 15081 1509 ) states that among the set of L2-normalized Dirichlet (or Neumann) 
Laplacian eigenfunctions, there is a sequence Uj,^ of density 1 (i.e., lim jk/k = 1), 

k^oo 

such that for any open subset V C 1^, one has 



lim 

k^oo 

V 



(this version of the theorem was formulated in [llOj. Marklof and Rudnick extended 
this theorem to rational polygons (i.e., simple plane polygons whose interior is con- 
nected and simply connected and all the vertex angles are rational multiplies of tt) 
[335j . Loosely speaking, {uj^^} is a sequence of non- localized eigenfunctions which 
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become more and more uniformly distributed over the domain (see jllQ[ 12011 1254] 
for further discussion and references). At the same time, this theorem does not pre- 
vent the existence of locahzed eigenfunctions. How large the excluded subsequence of 
(localized) eigenfunctions may be? In the special case of arithmetic hyperbolic man- 
ifolds, Rudnick and Sarnak proved that there is no such excluded subsequence [432 . 
This statement is known as the quantum unique ergodicity (QUE). The validity of 
this statement for other dynamical systems (in particular, ergodic billiards) remains 
under investigation [47l I159[ 1230] . The related notion of weak quantum ergodicity 
was discussed by Kaplan and Heller [268 . A classification of eigenstates to regular 
and irregular ones was thoroughly discussed (see [409, 495^ and references therein). 

There were numerous studies of Laplacian eigenfunctions in chaotic domains such 
as, e.g., Bunimovich stadium [69, 98i Ml ilMl UM EZSl [4831 HM., Sinai's bilhard 
[463, 464 . mushroom billiard [48^ .100] or hyperbolic billiard [2], illustrated on Fig. 



7.16 The literature on quantum billiards is vast, and we only mention selected works 
on the spatial structure of high-frequency eigenfunctions. McDonald and Kaufman 
studied the Bunimovich stadium billiard and reported a random structure of nodal 
lines of eigenfunctions and Wigner-type distribution for eigenvalue spacings |348[l349j . 
Bohigas and co-workers analyzed eigenvalue spacings for the Sinai's billiard and also 
obtained the Wigner-type distribution [81^. It means that eigenvalue spacings for 
these chaotic billiards obey the same distribution as that for random matrices from 
the Gaussian Orthogonal Ensemble. This is in a sharp contrast to regular billiards 
for which eigenvalue spacings generally follow a Poisson distribution. The problem of 
circular-sector and related billiards was studied (e.g., see |324j ). 

Polygon billiards have attracted a lot of attention, especially the class of rational 
polygons for which all the vertex angles are rational multiples of tt |325[ I326[ 1427] . 
As the dynamics in rational polygons is neither integrable nor ergodic (except several 
classical integrable cases such as rectangles, equilateral triangle, right triangles with 
an acute vertex angle 7r/3 or 7r/4), they are often called pseudo-integrable systems. 
Bellomo and Uzer studied scarring states in a pseudo-integrable triangular billiard 
and detected scars in regions which contain no periodic orbits [57]. Amar et al. gave 
a complete characterization of the polygons for which a Dirichlet eigenfunction can 
be found in terms of a finite superposition of plane waves ^ (see also |310j for 
experimental study). Biswas and Jain investigated in detail the 7r/3-rhombus billiard 
which presents an example of the simplest pseudo-integrable system [72] . Hassell et al. 
proved for an arbitrary polygonal billiard that eigenfunction mass cannot concentrate 
away from the vertices |229j (see also [336]). The level spacing properties of rational 
and irrational polygons were studied numerically by Shimizu and Shudo j457j . They 
also analyzed the structure of the related eigenfunctions [458 . 

Backer and co-workers analyzed the number of bouncing ball modes in a class of 
two-dimensional quantized billiards with two parallel walls [35 . Bunimovich intro- 
duced a family of simple billiards (called "mushrooms" ) that demonstrate a continuous 
transition from a completely chaotic system (stadium) to a completely integrable one 
(circle) [100 . Barnett and Betcke reported the first large-scale statistical study of 
very high-frequency eigenfunctions in these billiards [48 . Using nonstandard numer- 
ical techniques [47], Barnett also studied the rate of equidistribution for a uniformly 
hyperbolic, Sinai-type, planar Euclidean billiard with Dirichlet boundary condition, 
as illustrated on Fig. |7.17| This study brought a strong numerical evidence for the 
QUE in this system. The spatial structure of high-frequency eigenfunctions shown 
on Fig. 7.17 looks somewhat random. This observation goes back to Berry who con- 
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jectured that high-frequency eigenfunctions in domains with ergodic flow should look 
locally like a random superposition of plane waves with a fixed wavenumber [62 . This 
analogy is illustrated on Fig. 7.18 by Barnett |47]. O'Connor and co-workers analyzed 
the random pattern of ridges in a random superposition of plane waves j374j . 

Pseudo-integrable barrier billiards were intensively studied in a series of theoret- 
ical, numerical and experimental works by Bogomolny et al. [79l [80 . They reported 
on the emergence of scarring eigenstates which are related with families of classical 
periodic orbits that do not disappear at large quantum numbers in contrast to the 
case of chaotic systems. These so-called superscars were observed experimentally in 
a fiat microwave billiard with a barrier inside [80 . Wiersig performed an extensive 
numerical study of nearest-neighbor spacing distributions, next-to- nearest spacing 
distributions, number variances, spectral form factors, and the level dynamics [503 . 
Dietz and co-workers analyzed the number of nodal domains in a barrier billiard [156 . 

Tomsovic and Heller verified a remarkable accuracy of the semi-classical approx- 
imation that relates the classical and quantum billiards |483[ I484j . In some cases, 
eigenfunctions can therefore be constructed by purely semiclassical calculations. Li et 
al. studied the spatial distribution of eigenstates of a rippled billiard with sinusoidal 
walls |322j . For this type of ripple billiards, a Hamiltonian matrix can be found ex- 
actly in terms of elementary functions that greatly improves computation efficiency. 
They found both localized and extended eigenfunctions, as well as peculiar hexagon 
and circle-like pattern formations. Frahm and Shepelyansky considered almost circu- 
lar billiards with a rough boundary which was realized as a random curve with some 
finite correlation length. On a first glance it may seem that such a rough boundary in 
a circular billiard would destroy the conservation of angular momentum and lead to 
ergodic eigenstates and the level statistics predicted by random matrix theory. They 
showed, however, that there is a region of roughness in which the classical dynamics is 
chaotic but the eigenstates are localized |190j . Bogomolny et al. presented the exact 
computation of the nearest-neighbor spacing distribution for a rectangular billiard 
with a pointlike scatterer inside [78] . 

Prosen computed numerically very high-lying energy spectra for a generic chaotic 
3D quantum billiard (a smooth deformation of a unit sphere) and analyzed Weyl's 
asymptotic formula and the nearest neighbor level spacing distribution. He found 
significant deviations from the Gaussian Orthogonal Ensemble statistics that were 
explained in terms of localization of eigenfunctions onto lower dimensional classically 
invariant manifolds |410j . He also found that the majority of eigenstates were more 
or less uniformly extended over the entire energy surface, except for a fraction of 
strongly localized scarred eigenstates [411 . An extensive study of 3D Sinai's billiard 
was reported by Primack and Smilansky [408 . Deviations from a semi-classical de- 
scription were discussed by Tanner |481j . Casati and co-workers investigated how the 
interplay between quantum localization and the rich structure of the classical phase 
space infiuences the quantum dynamics, with applications to hydrogen atoms under 
microwave fields jll9H122] (see also references therein). 

A large number of physical experiments were performed with classical and quan- 
tum billiards. For instance, Graf and co-workers measured more than thousand first 
eigenmodes in a quasi two-dimensional superconducting microwave stadium billiard 
with chaotic dynamics [210 . Sridhar and co-workers performed a series of experi- 
ments in microwave cavities in the shape of Sinai's billiard |469[ 1470] . In particular, 
they observed bouncing ball modes and modes with quasi-rectangular or quasi-circular 
symmetry which are associated with nonisolated periodic orbits (which avoid the cen- 
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Fig. 7.17. Illustration of spatial distribution of the Dirichlet Laplacian eigenfunction |wmp 
(shown as density plots: larger values are darker) with m = 1, 10, 100, 1000 and m ~ 50000 jf^?/ (by 
A. Barnett, with permission) . 



Fig. 7.18. (Left) Density plot of a Dirichlet- Laplacian eigenfunction |wmP for m ~ 50000 
with the eigenvalue Am ~ 10^. There are about 225 wavelengths across the diagonal; (Right) 
Density plot of one sample from the ensemble of random plane waves with the same wavenumber 
and mean intensity, shown in a square region of space (with no boundary conditions) /^7/ (by A. 
Barnett, with permission). 



tral disk). Some scarring eigenstates, which are associated with isolated periodic 



orbits (which hit the central disk, see Fig. 7.16 3), were also observed. Kudrolli et al. 
investigated the signatures of classical chaos and the role of periodic orbits in the eigen- 
value spectra of two-dimensional billiards through experiments in microwave cavities 
5297' "298". The eigenvalue spectra were analyzed by using the nearest neighbor spacing 
distribution for short-range correlations and the spectral rigidity for longer-range cor- 
relations. The density correlation function was used for studying the spatial structure 
of eigenstates. The role of disorder was also investigated. Chinnery and Humphrey 
visualized experimentally acoustic resonances within a stadium-shaped cavity [131 . 
Bittner et al. performed double-slit experiments with regular and chaotic microwave 
billiards [73 . Chaotic resonators were also employed for getting specific properties of 
lasers (e.g., high-power directional emission or "Fresnel filtering") j2Q5[ I425 \ 

8. Other points and concluding remarks. This review was focused on the 
geometrical properties of Laplacian eigenfunctions in Euclidean domains. We started 
from the basic properties of the Laplace operator and explicit representations of its 
eigenfunctions in simple domains. After that, the properties of eigenvalues and their 
relation to the shape of a domain were briefly summarized, including Weyl's asymp- 
totic behavior, isoperimetric inequalities, and Kac's inverse spectral problem. The 
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structure of nodal domains and various estimates for the norms of eigenfunctions 
were then presented. The main Section [7| was devoted to the spatial structure of 
eigenfunctions, with a special emphasis on their localization in small subsets of a do- 
main. One of the major difficulties in the study of localization is that localization 
is a property of an individual eigenfunction. For the same domain, two consecutive 
eigenfunctions with very close eigenvalues may have drastically different geometrical 
structures (e.g., one is localized and the other is extended). One needs therefore fine 
analytical tools which would differently operate with localized and non-localized eigen- 
functions. In the review, we distinguished two types of localization, for low-frequency 
and high-frequency eigenfunctions. 

In the former case, an eigenfunction remains localized in a subset because of a 
geometric constraint that prohibits its extension to other parts of the domain. A 
standard example is a dumbbell (two domains connected by a narrow channel), for 
which an eigenfunction may be localized in one domain if its typical wavelength is 
larger than the width of the channel (meaning that an eigenfunction cannot "squeeze" 
through the channel). Such kind of "expulsion" from a channel is quite generic, as 
the analysis is applicable to domains with branches of variable cross-sectional profiles. 
It is important to note that a geometric constraint does not need to be strong (e.g., 
two domains may be separated by a cloud of point-like obstacles of zero measure). 
Another example is an elongated triangle, in which there is no "obstacles" at all. Low- 
frequency localization was found numerically in many irregularly-shaped domains, for 
both Dirichlet and Neumann boundary conditions. From a practical point of view, the 
low-frequency localization may find various applications, e.g., it is important for the 
theory of quantum, optical and acoustical waveguides and microelectronic devices, as 
well as for analysis and engineering of highly refiecting or absorbing materials (noise 
protective barriers, anti-radar coatings, etc). 

The high-frequency localization manifests in quantum billiards when a sequence 
of eigenfunctions tends to concentrate onto some orbits of the associated classical 
billiard. In this regime, the asymptotic properties of eigenvalues and eigenfunctions 
are strongly related to the underlying classical dynamics (e.g., regular, integrable or 
chaotic). For instance, the ergodic character of the classical system is refiected in the 
spatial structure of eigenfunctions. Working on simple domains, we illustrated several 
kinds of localized eigenfunctions which emerge for a large class of domains. We also 
provided examples of rectangular domains without localization. Although a number 
of rigorous and numerical results were obtained (e.g., quantum ergodicity theorem 
for ergodic billiards), many questions about the spatial structure of high-frequency 
eigenfunctions remain open, even for very simple domains (e.g., a square). 

Although the review counts more than five hundred citations, it is far from being 
complete. As already mentioned, we focused on the Laplace operator in bounded 
Euclidean domains and mostly omitted technical details, in order to keep the review 
at a level accessible to scientists from various fields. Many other issues had to be 
omitted: 

(i) Many important results for Laplacians on Riemannian manifolds or weighted 
graphs could not be included. In addition, we did not discuss the spectral properties of 
domains with "holes" [130l[l8a[287l[30l[3^ well as 

their consequences for diffusion in domains with static traps [21JJ I275( \276\ I369( 1485) . 
The behavior of the eigenvalues and eigenfunctions under deformations of a domain 
was partly considered in Sec. 6.1 and Sec. 7.5, while many significant results were 
not included (see [9ll [270l [2791 |423l |46l] and references therein). 



47 



(ii) There are important developments of numerical techniques for computing the 
Laplacian eigenbasis. In fact, standard finite difference or finite element methods rely 
on a regular or adapted discretization of a domain that reduces the continuous eigen- 
value problem to a finite set of linear equations ^ffl [El MB EH r299 , 435 . 
Since finding the eigenbasis of the resulting matrix is still an expensive computational 
task, various hints and tricks are often implemented. For instance, for planar polygo- 
nal domains, one can exploit the behavior of eigenfunctions at corners through radial 
basis functions in polar coordinates and the integration of related Fourier-Bessel func- 
tions on subdomains jl53[ I162[ I398j . Another "trick" is conformal mapping of planar 
polygonal domains onto the unit disk, for which the modified eigenvalue problem can 
be efficiently solved [37l[38]. Yet another approach, known as the method of particular 
solutions was suggested by Fox and co-workers [189" and later progressively improved 
[34| |49[ [67] . The main idea is to consider various solutions of the eigenvalue equa- 
tion for a given value of A and to vary A until a linear combination of such solutions 
would satisfy the boundary condition at a number of sample points along the bound- 
ary. One can also mention a stochastic method by Lejay and Maire for computing 
the principal eigenvalue |318) . The eigenvalue problem can also be reformulated in 
terms of boundary integral equations that reduces the dimensionality and allows for 
rapid computation of eigenvalues j332j . Kaufman and co-workers proposed a simple 
expansion method in which wave functions inside a two-dimensional quantum billiard 
are expressed in terms of an expansion of a complete set of orthonormal functions 
defined in a surrounding rectangle for which Dirichlet boundary conditions apply, 
while approximating the billiard boundary by a potential energy step of a sufficiently 
large size [273]. Vergini and Saraceno proposed a scaling method for computing high- 
frequency eigenmodes [496]. This method was later improved by Barnett and Hassell 
[50] (this reference is also a good review of numerical methods for high-frequency 
Dirichlet Laplacian eigenvalues). 

(iii) We did not discuss various applications of Laplacian eigenfunctions which 
nowadays range from pure and applied mathematics to physics, chemistry, biology 
and computer sciences. One can mention manifold parameterizations by Laplacian 
eigenfunctions and heat kernels [265 , the use of Laplacian spectra as a diagnostic 
tool for network structure and dynamics [350] , efficient image recognition and analysis 
[424", 426' "436' '437' , shape optimization and spectral partition problems [6^, 'UT, "1131 
[114 , 397 , 468 , 494 , computation and analysis of diffusion- weighted NMR signals [213^ 
I215j . localization in heterogeneous materials (e.g., photonic crystals) and the related 
optimization problem [EHl [T8THI841 [2621 [2631 [2851 l286l [2941 [2961 [438] , etc. 
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